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ABSTRACT 


An  integrated  circuit  tapped  delay  line  and  its 
applications  to  sampled  analog  processing  are  studied. 
The  theoretical  effort  includes  the  design  of  sampled 
analog  transversal  filters  using  the  techniques  developed 
for  digital  non-recursive  filters.   The  experimental  study 
includes  the  evaluation  of  the  new  Reticon  TAD- 12  tapped 
delay  line  and  its  applications  as  prewhitening ,  dewhiten- 
ing  and  bandpass  filters.   The  frequency  domain  response 
and  the  time  domain  impulse  response  are  evaluated.   The 
limitations  in  the  device  performance  and  the  deviations 
of  filter  performance  from  theoretical  analysis  are  investi 
gated. 
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I.   INTRODUCTION 

Signal  processing  covers  a  wide  range  of  application  in 
such  diverse  fields  as  biomedical  engineering,  acoustics, 
sonar,  radar,  speech  communication,  data  communication,  and 
many  others.   As  an  example,  one  may  wish  to  remove  interfer- 
ence, such  as  noise,  from  the  signal  or  to  modify  the  signal 
to  present  it  in  a  form  which  is  more  easily  interpreted  by 
an  expert.   As  another  example,  a  signal  transmitted  over  a 
communications  channel  is  generally  perturbed  in  a  variety  of 
ways,  including  channel  distortion,  fading,  and  the  insertion 
of  background  noise.   One  of  the  objectives  at  the  receiver 
is  to  compensate  for  these  perturbances .   In  each  case, 
processing  of  the  signal  is  required. 

Signal  proce-ssing  was  carried  out  typically  by  using 
analog  equipment  before  the  advent  of  digital  computers  which 
promoted  the  intensive  and  successful  development  of  digital 
signal  processing  in  the  past  decade.   At  present,  much  of 
the  digital  signal  processing  is  performed  on  general  purpose 
digital  computers  using  software.   However,  tremendous  progress 
has  also  been  made  in  other  types  of  signal  processors  using 
electronic  firmware  and  hardware.   Particularly,  using  elec- 
tronic technologies,  a  new  class  of  sampled  analog  processors 
is  being  developed  using  charge  transport  devices.   Their 
signals  are  analog  but  the  independent  variables  are  discrete 
which  can  be  either  time  variable,  spatial  variables  or  other 
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transform  variables.   The  signal  processing  implementations 
are  illustrated  in  Table  I.I. 

Signal  processors  can  be  grouped  into  two  classes: 
filtering  and  spectral  analysis.   This  classification  is 
highlighted  in  Figure  1.1.   All  of  these  signal  processing 
systems  require  the  basic  essential  functions,  such  as  cor- 
relation, convolution  and  transformation  (Fast  Fourier  trans- 
form). In  other  words,  they  are  shown  to  be  centered  around 
four  basic  mathematical  operations:   delay  (or  shift),  mul- 
tiplication, summation  and  generation  of  special  functions. 
These  operations  have  been  performed  by  digital  techniques 
and  it  was  probably  the  most  convenient  method.   However, 
the  new  sampled  analog  signal  processing  systems  only  require 
a  storage  for  discrete  time  analog  samples  and  not  discrete 
or   quantized  signal  amplitudes;  in  fact,  because  of  the 
computation  time,  frequency  selectivity,  and  part  cbunts,  it 
is  undesirable  to  quantize   the  signal.   Therefore,  discrete 
time  analog  devices  are  being  developed  to  perform  the  convo- 
lution, correlation,  etc. ,  without  the  need  to  convert  the 
signal  to  digital  format. 

There  are  several  distinct  advantages  of  using  the  discrete 
time  analog  devices  over  the  digital  techniques  to  perform  the 
same  algorithm: 

(1)  In  general,  the  discrete  time  device  system  can  perform 
the  same  computation  with  much  greater  speed  over  the 
digital  technique. 

(2)  The  device  stores  discrete  time  analog  amplitudes; 
therefore,  it  does  not  have  the  quantization  errors 
inherent  in  the  digital  system. 
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(3)   The  actual  multiplying  and  summing  operations  are 
easily  implemented,  and  with  fewer  components  than 
the  digital  system. 

To  illustrate  the  advantages  that  these  discrete  time 
devices  have  over  the  digital  technique,  a  simple  system 
shown  in  Figure  1-2  is  used.   This  is  a  model  of  a  simple 
single  pole,  first  order  discrete  time  filter.   To  implement 
this  system  with  a  discrete  time  device  such  as  the  SAD  (Serial 
Analog  Delay) ,  all  that  is  required  is  a  straight  substitution 
of  the  analog  device  for  the  delay  block,  an  operational 
amplifier  for  the  adder  and  a  potentiometer  proportional  to 
the  weighting  coefficient.   This  system  is  shown  in  Figure 
1-3. 

Comparing  this  system  implementation  to  the  digital  form 
in  Figure  1.4  demonstrates  the  difference  between  the  two 
systems.   It  requires  at  least  two  memory  registers,  one  to 
store  the  sampled  data,  the  other  to  store  the  multiplying 
constant.   In  the  figure  it  is  seen  that  more  devices  are 
required,  such  as  the  adder  and  the  costly  analog  to  digital 
converters,  as  well  as  the  digital  to  analog  converter  to 
complete  this  system. 

The  integrated  circuit  delay  devices  can  perform  in 
principle  most  of  the  necessary  algorithms  used  in  most  equip- 
ment for  signal  processing  today.   However,  due  to  the  primi- 
tive nature  of  these  new  classes  of  devices,  there  are 
limitations  in  some  fields  of  application. 

The  purpose  of  this  study  is  to  perform  the  evaluation  of 
Reticon  TAD-12  tapped  delay  lines  available  commercially  and 
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to  find  out  their  performances  and  limitations  in  order  to 
use  them  properly. 

TABLE  I . I   Signal  Processing  Implementation 


Signal 


Equipment 


Implementation 


Analog 


Digital  ( 


( ] 


Sampled 
Analog 


(    General  Purpose  Computer   Software 
Special  Purpose  Computer 
Mini/Microcomputer         Firmware 
(Microprocessor) 


Signal  Processing  Function 
Modules : 
Multiplier 
Correlator 


Convolver  (matched 
filter) 
Filter 

Discrete  Fourier 
Transformer 

Frequency  synthesizer 

etc. 


Hardware 
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Signal  Processors 


Filter 


Recursive 


Non- Recursive 


Digital 


Analog 


Correlator 


Matched  Filter 


Spectral  Analyzer 


Sampled  Analog 


Transversal  Filter 
Tapped  Delay  Line 


Reticon 
Device 


Prewhitening 
Dewhitening 


FIGURE  1.1   Classification  of  Signal  Processors 
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FIGURE  1.2  Equivalent  Circuit  of  First  Order  Discrete 
Time  Filter 
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FIGURE  1.3  Equivalent  Circuit  Implementation 
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FIGURE  1.4  Digital  Equivalent  Circuit  Implementation 
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II.   SAMPLED  ANALOG  NON- RECURSIVE  FILTER 

A.   PRINCIPLE 

The  transfer  function  of  a  digital  or  sampled  analog 
filter  can  be  expressed,  in  Z  domain,  by  one  of  the  follow- 
ing equations. 


Recursive  filter:   H(Z)  = 


ZN  a   Z"n 
L0   an  L 


,    „M  7   7-n 
1  m 


Non-Recursive  filter:   H(Z)  =  Eq  an  Z"n 


(2-1) 


(2-2) 


The  block  diagram  of  a  sampled  analog  recursive  filter 
is  illustrated  in  Figure  2-1: 


E.  (Z) 


Eout(Z> 


FIGURE  2.1   Schematics  of  a  Recursive  Filter 
(Tapped  Delay  Line) 
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-1-<g) 


out 


FIGURE  2.2     Schematic  of  a  Non-Recursive  Filter 
(Tapped  Delay  Line) 

In   the   case   of  causal   system  with   finite   duration  of 
impulse   response,    the   system   function   can  be  written   in   the 

form: 

N--1  _n 

(2-3) 


H(Z)    =nEQ   h(n)    Z 


{h(n)}  is  the  impulse  response  defined  over  the  time  interval 
0  <_  nT  <  (N  -  1)T,  H(Z)  is  a  polynomial  in  Z"1  of  degree  N  -  1 
Thus,  H(Z)  has  N-  1  zeroes  in  the  finite  Z-  plane. 

The  frequency  response  H(e^  )  is  obtained  by  substituting 

7  u    J  wT 
Z  by  eJ 
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uf    iwT.    N-l  ,  ,  .  -jwrtT  ro  .. 

H(eJ   J  =   E   h(n)  e  J  (2-4) 

n=0 

We  recall  that  any  finite  duration  sequence  is  completely 
specified  by  N  samples  of  its  Fourier  Transform  [1] ,  so  that 
the  design  of  a  Finite  Impulse  Response  Filter  may  be 
accomplished  by  finding  either  its  impulse  response  coeffi- 
cients or  N  samples  of  its  frequency  response.   Both  methods 
are  discussed  in  the  following  section. 

B.   LINEAR  PHASE  FINITE  IMPULSE  RESPONSE  FILTER 

In  many  applications,  e.g.,  speech  processing,  data 
transmission,  it  is  desirable  to  design  filters  to  have 
linear  phase.   In  this  way,  signals  in  the  passband  of  the 
filter  are  reproduced  exactly  at  the  filter  output  except 
for  a  delay  corresponding  to  the  slope  of  the  phase.   One 
of  the  most  important  features  of  FIR  systems  is  that  they 
can  be  designed  to  have  exactly  linear  phase.   The  unit 
impulse  response  for  a  causal  FIR  system  with  linear  phase 
has  the  property  that: 

h(n)  =  h(N  -  1  -  n)  (2-5) 

To  see  that  this  condition  implies  linear  phase,  we  write 
equation  2-4  as: 

(N/2)-l       _n     N-l 

H(2)=  I         h(n)  Z  n  +    *   h(n)  Z  n 
n=0  n=N/2 

(N/2)-l  (N/2)-l 

=   n^O   Mn)Z    +  nL     h(N-l-n)2-CN-l-n) 
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where  N  is  assumed  to  be  even.   Using  equation  2-5,  one  can 

write: 

CN/2)-l       _n    -fN-i-n-) 

H(Z)  =    Z     h(n)[Z  n  +  Z  lIN  -1  nj]  (2-6) 

n=0 
If  N  is  odd,  it  is  shown  that 

H(Z)  =  ^^^hdOIZ""  +  Z-CN-l-n)]  +h(Nii)z-[(N-l)/2] 
n=0  z 

(2-7) 

Evaluating  equations  2-6  and  2-7  for  Z  =  e-*w,  one  obtains, 
for  N  even, 

ur   jw.     -jw[(N-l)/2]  ,  ^N{2)"1  ~  f   .      r  r    N-1A11 
H(eJ  )  =  e  J  LV   y  J  {   E    2h(n)  cos  [w  (n  -  — s— )  J  } 

n=0  L 

(2-8) 

and  for  N  odd, 


H(e^)  =  e-WCN-l)/*]  {h(% 


[(N-3)/2]-l 
+     E        2h(n)  cos  [w (n  -  %&  ] }   (2-9) 
n=0  z 


In  both  cases  the  sums  in  brackets  are  real,  implying  a 
linear  phase  shift  corresponding  to  delay  of  (N-l)/2  samples. 
For  N  odd,  the  phase  shift  corresponds  to  an  integer  number 
of  samples  delay.   For  even  N,  the  delay  is  an  integer  plus 
one-half  of  sampling  period.   This  distinction  between  odd 
and  even  values  of  N  is  often  of  considerable  importance  in 
design  and  realization  of  FIR  filter.   Examples  of  impulse 
response  having  linear  phase  are  shown  in  Figure  2-3. 
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FIGURE  2.3     Typical  Impulse  Response  for  Linear  Phase  FIR  Filters 

C.      DESIGN  TECHNIQUE   OF    FIR  FILTERS 

1 .      Determination  of   Impulse   Response 

Since   the    frequency   response   of  H(e->    )    of   any   digital 
or   sampled   analog   filter   is   periodic    in   frequency,    it   can  be 
expanded   in   Fourier   series   of   the    form: 


H(eJW)    =      Z        h(n)    e 
n=-oo 


-jwn 


(2-10) 


where  the  sequence  h(n)  playing  the  role  of  the  "Fourier 
coefficients,"  i.e., 


jw.   jwn 


h(n)  =  j-   /"   H(eJW)  e 


dw 


(2-11) 
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One  can  easily  recognize  that  the  coefficients  of  the 
Fourier  series,  h(n) ,  are  identical  to  the  impulse  response 
of  a  digital  filter.   There  are  two  difficulties  with  the 
representation  of  (2-10)  for  designing  FIR  filters.   First, 
the  filter  impulse  response  is  infinite  in  duration  since  the 
summation  in  (2-10)  extends  to  ±°°  .   Second,  the  filter  is 
unrealizable  because  the  impulse  response  begins  at  -°°;  i.e., 
no  finite  amount  of  delay  can  make  the  impulse  response 
realizable. 

One  possible  way  of  obtaining  an  FIR  filter  that 
approximates  H(e^  )  would  be  to  truncate  the  infinite  Fourier 
series  (2-10),  at  n  =  ±  M.   But  direct  truncation  of  the 
series  leads  to  the  well  known  Gibbs  phenomenon.   However, 
the  non-uniform  convergence  phenomenon  can  be  moderated 
through  the  use  of  a  less  abrupt  truncation  of  the  Fourier 
series  by  applying  the  so-called  "windowing"  or  "weighting 
function. " 

2 .   Effect  of  Window 

A  more  successful  way  of  obtaining  an  FIR  filter  is 
to  use  a  finite  weighting  sequence  W(n) ,  called  window,  to 
modify  the  Fourier  coefficients  h(n)  in  (2-10)  to  control  the 
convergence  of  the  Fourier  series.   The  technique  of  windowing 
is  illustrated  in  Figure  2-4.   At  the  top  of  this  figure  is 
shown  the  desired  periodic  frequency  response  H(e-JW)  and  its 
Fourier  series  coefficients  (h(n)}.   The  next  row  shows  a 
finite  duration  weighting  sequence  W(n)  with  Fourier  trans- 
form W(eJ  ).   W(e-J  ),  for  most  reasonable  windows,  consists 
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of  a  central  lobe  which  contains  most  of  the  energy  of  the 
windows  and  side  lobes  which  generally  decay  rapidly.   To 
produce  an  FIR  approximation  to  H(e^  ) ,  the  sequence 

*   ,     h(n)  .  W(n)  ,  -M  <  n  <  M  (2-12) 

u  ,   elsewhere 

is  formed.   The  third  row  of  Figure  2.4  shows  h(n)  and  its 
Fourier  transform  H(e^w) ,  which  is  the  convolution  of 
H:(e^  )  and  W(e^  ),  since  h(n)  is  the  product  of  the  sequences 
h(n)  and  w(n) .   The  last  row  of  Figure  2-4  shows  the  realiz- 
able sequence  g(n),  which  is  a  shifted  version  of  h (n)  and 
may  be  used  as  the  desired  filter  impulse  response. 

As  seen  in  the  example  of  Figure  2-4,  there  are 
several  effects  on  the  resulting  frequency  response  of  win- 
dowing the  Fourier  coefficients  of  the  filter.   A  major  effect 
is  that  discontinuities  in  H(e^w)  become  a  transition  band 
between  values  on  either  side  of  the  discontinuity.   Since 
the  final  frequency  response  of  the  filter  is  the  circular 
convolution  of  the  ideal  frequency  response  with  the  window's 
frequency  response,  it  is  clear  that  the  width  of  these  tran- 
sition bands  depend  on  the  width  of  the  main  lobe  of  WCe-*  )  . 
A  secondary  effect  of  windowing  is  that  ripple  from  the  side 
lobes  of  W(e^  )  produces  approximation  errors  (ripple  in  the 
resulting  frequency  response)  for  all  w. 

There  are  many  windows  proposed  that  approximate  the 
desired  characteristics.   But,  in  general,  the  desirable 
window  should  have  the  following  characteristics: 
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FIGURE  2.4  Illustration  of  windowing 
H(e^w)  Desired  periodic  frequency  response 

h(n)   Its  Fourier  series  coefficients 
w(n)    Finite  duration  weighting  sequence 
w(eJW)  Its  Fourier  transform 
h(n)  =  h(n)  .  w(n) 
g(eJw)-  H(eJW)  *  W(ejW) 
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1.  Small  width  of  main  lobe  of  the  frequency  response  of 
the  window  containing  as  much  of  the  total  energy  as 
possible. 

2.  Side  lobes  of  the  frequency  response  that  decrease  in 
energy  rapidly  as  w  tends  to  it. 

In  the  design  of  the  filters  in  this  study,  the 

Hamming  window  is  chosen.   It  is  specified  by  the  equation: 

r  0.54  -  0.46  cos  (!-t),   0  <  n  <  N-l 
w(n)  =  {  N"1      -   - 

0  ,   elsewhere     f2-131 

Figure  2-5  illustrates  the  Hamming  window  and  the  rectangular 
window. 


W(n) 


Rectangular 


n 


FIGURE  2.5     Hamming  Window  and  Rectangular  Window 


3.      Frequency   Sampling   Technique 

Figure    2-6   shows    an   arbitrary   frequency   response 
(solid  curve)    that   one   wants    to   approximate   and  a   sequence   of 
N   frequency   samples   H    (Figure    2-7)    that   can  be   represented   as 


H(k)    =    |H(k) |e 


j6(k) 


k   =    0,    1.  .  . ,    N-l 
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FIGURE  2.7     Frequency  Samples 


Using   the    inverse   discrete   Fourier   transform,    a 
finite   duration   impulse   response   can  be   determined   from  H(k) 


N-1 


h(n)    .   1      S      H(k)    e^2^/N^nk 


N 


k=0 


(2-14) 


where  H(k)    can   also   be  written   in   the    form: 

H(k)    =   Y   h(n)    e-H2iT/N)nk 

n=0 


(2-15) 
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D.   APPLICATIONS 

1.   Bandpass  Filter  Design  Algorithm 

In  the  analog  or  continuous  time  system,  the  bandpass 
filter  has  the  transfer  function  H(f) . 


f   f   f 
xl  r0  r2 


FIGURE  2.8   Bandpass  filter  frequency  response 


H(f)  =  { 


A,     f±   <    f  <  f2 

0,     otherwise 


In  the  discrete  time  system,  a  rectangular  bandpass 
filter  has  the  frequency  response  HCe-'™)  as  shown  in  Figure 
2-9.   It  is  periodic  of  period  2tt  . 


w 


H(eJws 


A>    lwil  < 


wj|  <  |w2 


0,    otherwise 


w< 


with  f0  = 

L  2lT 


sampling  frequency. 
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FIGURE  2.9     Discrete  Time  Bandpass  Filter  Frequency  Response 

Basing  upon  the  periodic  nature   of  a   discrete   time 
filter   frequency  response,    the   impulse   response   can  be   com- 
puted by   finding   the   Fourier  coefficient,   using   equation   2-11 


h(n)    -   £  /    2   A  J™  dw 

1 
1         W2 

=  =■  /        A    (cos  wn   +   j    sin  wnj    dw. 


(2-16) 


Since  only  real  term  may  be  implemented,  the  series  coeffi 
cients  are: 

Case  I :  n  --h   0 

w. 


h(n)  =  \   /  2  A  cos  (wn)  dw 

=   A.   rsin   2irf      -    sin   2irf   ) 
-rrn    v  2  i 


(2-17) 


Case  2 
Let 


n  =  0 


w 


w 


=  2  it  (fo  +  &f) 

=  2  IT  (fo  +  ^f) 


One  can  write  using  equation  (2-16) 
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to: 


w 

h(n)    =  —  f        eJ        dw 

it      w 

=   JL.   [eJW2n    -    ejwin] 

=   J±-  ei2irf°    (sin   TrnAf) 

=    2A  Af51"  ^An| f    (cos    2TTnf0   +   j    sin   2tt  n  £0) 

(2-18) 
When  n  tends  to  0,  the  real  part  of  equation  (2-18)  tends 


h(0)  =  2AAf(l.)(l.)  =  2AA  f  (2-19) 


Window  weighting  functions  are  used  to  modify  the 
coefficients  in  order  to  improve  the  filter  characteristics. 
If  the  Hamming  window  is  used,  the  weighted  coefficients 
A(n)  become 

A(n)  =  h(n)(0.54  +  0.46  cos  (^tt))  (2-20) 

A  computer  program  shown  in  Appendix  C  was  established 

to  design  the  non-recursive  bandpass  filter.   It  was  based 

upon  the  equations  2-17,  2-19,  2-20. 

2 .   Prewhitening  and  Dewhitening  Filters  in 
Communication  System 

The  prewhitening  and  dewhitening  filters  are  also 
called  preemphasis  and  deemphasis  filters  in  the  communication 
field. 

Suppose  one  desires  to  transmit  a  baseband  signal 
using  FM  modulation  and  requires  the  best  possible  signal  to 
noise  ratio.   One  possible  solution  is  to  raise  the  level  of 
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modulating  baseband  signal  to  the  maximum  extent  possible 
in  order  to  modulate  the  carrier  as  vigorously  as  possible. 
But  the  modulating  signal  level  may  be  raised  only  until  the 
distortion  exceeds  a  specified  value  (maximum  allowable  value) 

However,  the  baseband  signal  happens  to  be  an  audio 
signal;  it  turns  out  that  something  further  can  be  done.   An 
audio  signal  usually  has  the  characteristic  that  its  power 
spectral  density  is  relatively  high  in  the  low  frequency 
range  and  falls  off  rapidly  at  higher  frequencies.   For  ex- 
ample, speech  has  little  power  spectral  density  above  about 
3  kHz.   Music  extends  farther  into  the  high  frequency  range, 
the  feature  still  persists  that  most  of  its  power  is  in  the 
low  frequency  region.   As  a  consequence,  when  one  examines 
the  spectrum  of  the  sidebands  associated  with  a  carrier  which 
is  frequency  modulated  by  an  audio  signal,  one  finds  that  the 
power  spectral  density  of  the  sidebands  is  greatest  near  the 
carrier  and  relatively  small  near  the  limits  of  the  allowable 
frequency  band  allocated  to  the  transmission.   The  manner  in 
which  one  may  take  advantage  of  these  spectral  features, 
which  are  characteristic  of  audio  signals,  in  order  to  improve 
the  performance  of  an  FM  system  is  shown  in  Figure  2-10. 

At  the  transmitting  end,  in  Figure  2-10,  the  baseband 
signal  m(t)  is  not  applied  directly  to  the  FM  modulator  but 
is  first  passed  through  a  filter  of  transfer  characteristic 
Hp (w) ,  so  that  the  modulating  signal  is  mp (t) .   The  modifi- 
cation introduced  into  the  baseband  signal  by  the  first 
filter  is  undone  by  the  receiver  filter  which  follows  the 
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FIGURE  2.10  Preemphasis  and  Deemphasi3  in  an  FM  system 
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discriminator  and  has  transfer  characteristic  1/H  (w)  .   The 
noise  passes  through  only  the  receiver  filter  which  may  be 
used  to  suppress  the  noise  to  some  extent. 

The  selection  of  the  transfer  characteristic  Hp(w) 
is  based  on  the  following  considerations.   At  the  output  of 
the  demodulator,  the  spectral  density  of  the  noise  increases 
with  the  square  of  the  frequency,  as  shown  in  Figure  2-11. 
It  is  given  by  the  equation: 


Gn(f)  = 


°Q  w2 


i! 


(2-21) 


where  a  =  constant  proportional  to  limited  amplitude 

of  the  carrier  obtained  with  the  hard  limiter 
in  the  demodulator 

A  =  amplitude  of  the  carrier 

n  =  power  spectral  density  of  the  white  noise 
over  the  range    |f|  <  y  • 


|on(f)l 


FIGURE  2.11   Power  Spectral  Density  at  the  Output 
of  an  PM  Demodulator 
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Hence  the  receiver  filter  will  be  most  effective  in 
suppressing  noise  if  the  response  of  the  filter  falls  off 
with  increasing  frequency.   In  such  a  case  the  transmitter 
filter  must  exhibit  a  rising  response  with  increasing  fre- 
quency.  The  premodulation  filtering  in  the  transmitter,  to 
raise  the  power  spectral  density  of  the  baseband  signal  in 
its  upper  frequency  range,  is  called  preemphasis.   The  fil- 
tering at  the  receiver  to  undo  the  signal  preemphasis  and 
to  suppress  noise  is  called  deemphasis. 

Referring  to  Figure  2-10,  one  requires  that  the  nor- 
malized power  of  the  baseband  signal  m(t)  must  be  the  same 
as  the  normalized  power  of  the  preemphasized  signal  m  (t) . 
If  G  (f)  is  the  power  spectral  density  of  m(t) ,  the  density 

of  m  ft)  is  |H  (f ) I 2  G  (f) ,  and  one  requires  that: 
p       '  p    '    m.  n 

fm  fm 

Pm ■  -  /   G  (f)  df  =  /   |H  (f)  |2G  (f)  df  (2-22) 

m    .£m  m^        -fml  p^  J  '  m"-  J 

where 

fm  =  maximum  frequency  of  modulating  signal. 

In  the  absence  of  deemphasis,  output  noise  is: 


8tt2      q2n 
T--aT- 


No    =    £Ji_  J±Jl    £m: 


With   deemphasis,    the   output   noise    is: 
2  fm  -, 

Nod  -  <£o   ^  '-a.*1  ifyni2d£  v-iv 
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No 
The   ratio,    R  =   « —   ,    is 

od 

R  =      t=      fm3/3 (2-24) 

/*m   f2df/|Hp(f)|2 

Since  the  signal  itself  is  unaffected  in  the  overall 
process,  the  quantity  R  is  the  ratio  by  which  preemphasis- 
deemphasis  improves  the  signal  to  noise  ratio. 

The  prewhitening  filter  has  the  transfer  character- 
istic illustrated  in  Figure  2-12.   The  dewhitening  filter 
has  the  transfer  characteristic  illustrated  in  Figure  2-13. 
The  relation  between  these  two  filters  is: 

where 

H  (e-'  )   =  transfer  function  of  prewhitening  filter 

i  w 
Hn(e   )   =  transfer  function  of  dewhitening  filter. 

3 .   Prewhitening  Filter  Design  Algorithm 

a.   In  the  analog  or  continuous  time  system,  the 
prewhitening  filter  has  the  transfer  characteristic  shown 
in  Figure  2-12  and  has  the  transfer  function, 


R  ,    £  c  f  0 

H  (f)  = 

R  +  A(f  -  f0)  ,    f0  <   f|  <  f 


(2-26) 
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f   f  normalized 
P 


FIGURE  2.12  Prewhitening  Filter  Transfer  Function 


yr)t 


R 


f^  f     f  normalized 

0  p 

FIGURE  2.13     Dewhitening  Filter  Transfer  Function 


b.      In   the   sampled  analog   or  digital    system,    the 
frequency  response   of   the  prewhitening   filter    is   periodic 
It    is    illustrated   in   Figure    2-14. 
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(ej2W)l 


normalized  frequency 
FIGURE  2.14  Amplitude  Frequency  Response  of  Sampled 

Analog  Prewhitening  Filter 


The   frequency   F    is   normalized  with   respect   to 
the   sampling    frequency   F    . 

The   design  of  a   linear  phase  prewhitening   filter 
consists   of   the   calculation  of   the    impulse   response  which 
can  be   computed  using    (2-11) 


h(n)    =  j-  f     H(eJ    )    eJ 


dw 


which  becomes,  using 


w 

Wc 


=    F 


F  Fn 

h(n)    =    2    /   °    R  eJ2™F   dF    +    2/c       [R   +   A(F-Fj]eJ2TrnFdF 
O  r  O 

O 


h(n)  = 


U 


(2-27) 


Since  h(n)  is  real,  one  can  compute  only  the  real 
part  of  U  and  V. 
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Fs 
Case   1:      n  f   0    ,  F     =  -j- 

Fs  R 

Re[U]    =    2R  /        cos    (27rnF)dF   =   —  sin   27rn  FQ 
i    j  0  ^  y  irn  ° 

F  F 

Re[V]    =    2R  /  p   cos    (27rnF)dF   +    2A  /  P  (F  -  F0)  cos  (2TrnF)dF 
Fo  Fo 

=  — -   [sin    (2irn  F_.)    -    sin    (2irn   Fn)] 
irn    L  v  p  °  ° 

+    2A    [4   2n2     (cos    2irnFp    -    cos    2iTnF0) 


♦2 


AF 


FD 

—  sin   27rnFn    -   -^L—  sin   27rnFn] 
n  P        2irn  UJ 


+   — —   [sin   2irnFn    -    sin   2irnFn] 
irn      L  u  P 


h(n)    =   Re    [U]    +   Re    [V]    =    7   ft  2    [cos    27rn   Fp  -  cos    2irn  F0  ] 

(2-28) 
Case    2:      n   =    0    ,  Fp   =    Fs/2 

With  n  =  0,  (2-27)  becomes: 

Fp  Fp 

h(o)    =    2    /        RdF  +  2  /      [R  +   A(F  -  FQ)  ]  dF 

Fo 

h(0)   =  ~TPj  [R  +  2F^  (Fp  "  F°)2]  (2'29) 

Taking  into  account  the  constraint  of  the  linear 
phase  FIR  filter: 

h(n)  =  h  (N  -  1  -  n) 
and  using  (2-28)  and  (2-29),  a  computer  program 
is  prepared  to  calculate  the  impulse  response  of 
the  filter,  a  Hamming  window  function  is  used  to 
improve  the  amplitude  and  phase  response  of  the 
prewhitening  filter. 
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4.   Dewhitening  Filter  Design  Algorithm 

In  the  design  of  the  complementary  dewhitening  filter, 
one  can  use  the  frequency  sampling  technique.   Using  (2-15) 
and  the  given  impulse  response  of  the  prewhitening  filter, 
the  frequency  samples 

{H  (k)}   ,    k  =  0,  1 ,  N-l 

can  be  computed  and  stored  in  the  memory  of  the  computer. 
Then  the  frequency  samples  of  dewhitening  filter 
H]}(k)  can  be  computed  as: 

HD(k)  =  H~7kT   •    k  =  0,  1,....,  N-l 

Using  the  inverse  discrete  Fourier  transform  given 
by  (2-14)  ,  the  finite  impulse  response  of  the  dewhitening 
filter  can  be  obtained. 

Using  an  appropriate  window,  such  as  Hamming  window, 
the  frequency  response  of  the  dewhitening  filter  can  be 
improved  also. 


44 


III.   EXPERIMENTAL  EVALUATION  OF  TAD- 12 


A.   SYSTEM  CONFIGURATION 

Reticon  Corporation  markets  a  Tapped  Analog  Delay  Device 
which  provides  12  output  taps  with  the  following  time  delay 
configuration: 


In 


T 


2T 


2T 


2T 


2T 


12 


) 


FIGURE  3.1  Iqulralent  Circuit  of  TAD-12 


The  first  delay  section  is  delayed  by  one  clock  period, 
the  rest  have  delay  time  of  two  clock  periods. 

The  delay  lines  are  structurally  organized  as  shown  in 
Figure  3-2.   One  delay  line  contains  N  capacitive  storage 
elements  with  two  I/O  access  switches.   One  access  switch 
connects  the  storage  element  to  a  common  input  line  and  the 
second  switch  connects  it  to  the  common  output  line.   One 
register  and  its  associated  multiplexing  switch  sequentially 
time  samples  the  analog  input  signal  onto  the  "N"  capacitive 
storage  elements.   Information  is  read  onto  the  nth  element 
while  simultaneously  reading  out  the  n   +1  element. 
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FIGURE  3.2  Organization  of  a  Single  Delay  Line 
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Each  element  has  an  output  buffer  amplifier  to  provide 
a  zero  order  hold  of  the  sampled  signal.   This  produces  a 
sampled  and  held  output. 

The  test  circuit  for  TAD-12  is  shown  in  Figure  3-3.   It 
consists  of  a  pulse  shaper  circuitry  and  a  DMOS  flip-flop. 
Because  the  tap  outputs  are  at  +5  volts  bias  level,  due  to 
the  TAD-12  circuitry,  the  resistor  2  kilo-ohms  of  the  oper- 
ational amplifier  non- inverting  input  is  connected  to  a  +5 
volts  rather  than  ground.   The  output  of  the  summer,  there- 
fore, rides  on  a  5  volts  bias  level. 

B.   FREQUENCY  DOMAIN  EVALUATION 

Using  the  circuit  board  shown  in  Figure  3-3,  the  frequency 
dependance,  the  non-uniformity,  the  harmonic  distortion  are 
evaluated  at  different  sampling  frequencies  for  three  differ- 
ent tap  resistances: 

Rk  ■   10  Kfl   ,         100   Kfl    ,        1  MJ2    . 

1 .   Frequency  Roll-off 

The  frequency  characteristic  of  the  sampled  and  hold 

s  in  x  7T  f 

process  is  of  the  shape  with  x  =  -z—  . 

x  rs 

At  sampling  frequency  f  =  500  KHz,  for  R.    =  10  Kfi, 

s  l.  ap 

the  roll-off  is  close  to  the  theory.   When  the  tapping 

resistances  R.    is  increased,  the  roll-off  is  more  serious 

tap  ' 

and  wriggles  around. 

Figure  3-4  illustrates  the  magnitude  of  the  frequency 
response  of  Tap  #1,  Tap  #6  and  Tap  #11  outputs  using  10  kilo- 
ohms  resistor.   The  solid  curve  represents  the  theoretical 
response.  47 
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FIGURE  3.3  Test  Circuit  for  Ret  icon  TAD- 12  Evaluation 
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From  this  experiment  it  is  found  that  the  TAD- 12 

output  has  rather  complicated  frequency  variations  although 

s  in  x 
they  follow  the  theoretical  relation  roughly  if  the 

tapping  resistor  is  10  Kilo-ohms  and  the  sampling  frequency 
less  than  2.5  MHz.   The  loading  effect  is  much  more  pronounced 
if  the  tapping  resistors  are  equal  or  greater  than  100  Kilo- 
ohms. 

2.   Non-Uniformity  (Spatial  or  Fixed  Pattern  Noise) 

The  non-uniformity  of  different  taps  outputs  varies 
depending  on  the  tap  resistance,  R,  ,  the  sampling  frequency, 
f s ,  and  the  input  signal  frequency,  f.   The  general  result 
at  sampling  frequency  of  500  KHz  is  the  following: 

-  For  f  <  10  KHz,  the  non-uniformity  is  within  101 

-  For  f  >  10  KHz,  the  non-uniformity  can  deteriorate 
to  15%  -  50%  range. 

Figures  3-10  through  3-17  illustrate  the  non-uniformity 
result  of  sampling  frequency  of  500  KHz,  for  different  input 
signal  frequencies  and  for  different  tap  resistors: 

Rtap  =  X  Kf2,   10  K^'   10°  K^»   X  m' 

At  high  sampling  frequencies,  using  tap  resistances 
less  than  100  Kfi,  the  uniformity  of  tap  output  can  be  improved 
at  higher  input  signal  frequencies. 

At  fs  =  0.5  MHz,  R,  =  10  Kfi  and  input  signal  frequencies: 

-  f  <  50  KHz,  the  non-uniformity  is  within  10% 

-  f  >  50  KHz,  the  non-uniformity  varies  from  13%  to  25%. 
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At  fs  =  2.5  MHz,  R,  =  100  Kfi  and  input  signal  frequencies: 

-  f  <  80  KHz,  the  non- uniformity  is  within  101 

-  f  >  80  KHz,  the  non-uniformity  varies  from  11%  to  70% 

At  f  =2.5  MHz,  R,  =  1  Mil,  and  input  signal  frequencies: 

-  f  <  30  KHz,  the  non-uniformity  is  within  201 

-  f  >  30  KHz,  the  non-uniformity  is  greater  than  30%. 

Figures  3-18  through  3-28  illustrate  the  non-uniformity 
data  at  sampling  frequency  equal  to  2.5  MHz,  for  different 
typical  tap  resistors  and  for  different  input  signal  fre- 
quencies . 

Figures  3-24  through  3-40  illustrate  the  pictorial  pres- 
entation of  Reticon  TAD- 12  performance  at  different  sampling 
frequencies,  using  different  tap  resistances. 

3.   Temporal  Noise  Behavior  of  Different  Taps 

The  temporal  noise  of  different  taps  is  studied  at 
the  sampling  frequencies  higher  or  equal  to  2.5  MHz.   It 
varies,  depending  on  the  sampling  frequencies  and  the  tap 
resistances  used.   However,  it  should  be  pointed  out  that 
the  tap  #1  output  has  little  noise  at  all  sampling  frequen- 
cies and  tapping  resistors  used.   It  is  shown  in  Figures  3-24, 
3-25,  3-26,  3-28,  3-40.   The  general  temporal  noise  behavior 
of  other  taps  is  shown  in  Table  III. I  and  Figures  3-24,  3-25, 
3-27,  3-35,  3-38. 
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TABLE  III.  I  Noisiest  Tap  Outputs 

Rtap 
f  1  Meg-ohm    100  Kilo-ohms    10  Kilo-ohms 

2.5  MHz       Tap  #10       Tap  #11         Tap  #12 
3.5  MHz       Tap  #10       Tap  #4  Tap  #2 

The  temporal  noise  of  different  taps  is  not  serious 
at  low  sampling  rate  (less  or  equal  to  1.25  MHz).   Figure 
3-41  illustrates  the  output  of  tap  #1  and  #12  at  sampling 
frequency  at  fs  =  1.25  MHz  using  1  Meg-ohm  tapping  resistors. 
4.   Harmonic  Distortion 

The  second  harmonic  distortion,  data  is  measured  with 
the  input  frequency  10  KHz  for  different  input  amplitudes  at 
different  sampling  frequencies: 

fs  =  500  KHz,  1.25  MHz,  2.5  MHz 

and  using  different  tap  resistors: 

R<_    =  10  Kfl,  100  Kfl,  1  Mfi. 
tap  ' 

The  curves  for  the  quantity  20  log  y  fundamental 
are  plotted  and  shown  in  Figures  3-42,  3-43,  3-44.   The  general 
results  of  the  above  measurement  are  summarized  in  Table  3. II. 

The  loading  effect  on  the  dynamic  range  of  the  TAD-12 
is  very  pronounced  at  high  sampling  frequency  (2.5  MHz).   In 
order  to  get  the  second  harmonic  amplitude  less  than  1%   of 
the  fundamental,  one  has  to  use  less  than  .12  volt  rms  input 
voltage  if  using  1  Mft  tap  resistor  and  less  than  2  volts  rms 

input  voltage  if  using  10  Kfl  tap  resistor. 
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Number 1 1 
of  Taps' 


.04 


20%  variation 


.045 
Tap  Output  Amplitude 


T    Y 


-* »-» 


.05  volts 


FIGURE  3.10  Non  Uniformity  Data,  f    »  500  KHz 

samp 

R.   -  1  Kilo-ohm.  f.    .-  45  KHz 
tap  input 


Number 
of  Taps 
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l   I 


.03 


55%  variation 


.04 
Tap  Output  Amplitude 
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FIGURE  3.11  Non  Uniformity  Data,  f    »  500  KHz 

samp 

R.   »  1  Kilo-ohm.  f .    .-  60  KHz 
tap  input 


58 


Number 

i 

of  Taps 

25$  variation 

3 

• 

« 

i 

2 

• 

> 

3 

1 — 

< 

i 

< 

< 

I 

< 

< 

* 

—+ 

.1 


.11         .12 
Tap  Output  Amplitude 


.13  rolts 


FIGURE  3.12  Non  Uniformity  Data,  *saiip-  500  KHz 
Rtap-  10  K2lo-ohBls-   finPut"  5°  mz 
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FIGURE  3.13  Non  Uniformity  Data,   f    -  500  KHz 

98&H 

R^     =  10  Kilo-ohms,      f .        .-  100  KHz 
tap  input 
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22$  variation 
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Tap  Output  Amplitude 
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FIGURE  3.14  Non  Uniformity  Data,   f    -  500  KHz 
Htap"  10°  Kil°-ohBS-   f input"  5°  m 


Number 
of  Taps 


45$  variation 


.13  .H   .15   .16   .17   .18   .19  volts 
Tap  Output  Amplitude 


FIGURE  3.15  Non  Uniformity  Data,   f    -  500  KHz 

SCUHJU 

R.   -100  Kilo-ohms.   f.    .-  60  KHz 
tap  input 
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Tap  Output  Amplitude 
FIGURE  3.16  Non  Uniformity  Data,  f 


11.3$  variation 
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R.   -  1  Meg-ohm.  f .   . 
tap  input 
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samp 
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FIGURE  3.17  Non  Uniformity  Data,   f    -  500  KHz 

samp 

R,   »  1  Meg-ohm.  f.    .»  120  KHz 
tap      °       input 
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FIGURE  3.18  Non  Uniformity  Data,   f    ■  2.5  MHz 

samp 

Rx  -  10  Kilo-ohms,  f .    .-  50  KHz 
tap  input 


Number 
of  Taps 
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8  . 
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25$  variation 
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FIGURE  3.19  Non  Uniformity  Data,   f    =2.5  MHz. 

samp 

Rx  -  10  Kilo-ohms.   f.    .-  350  KHz 
tap  input 


.5  volts 


>  » 


,5  volts 
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10.8$  variation 
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FIGURE  3.20  Non  Uniformity  Data,  f    =2.5  MHz 

samp 

R.   »  100  Kilo-ohms,   f .    .-  80  KHz 
tap  input 


Number 
of  Taps 

5  ' 
4  ' 
3 


.1 


50$  variation 


.2         .3 
Tap  Output  Amplitude 

FIGURE  3.21  Non  Uniformity  Data,  faamp-  2.5  MHz 
R   -  100  Kilo-ohms.   finput-  300  KHz 


.4  volts 
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27.3%  variation 
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.08        .12 
Tap  Output  Amplitude 
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FIGURE  3.22  Non  Uniformity  Data,  f    -  2.5  MHz 

samp 

R.   =■  1  Meg-ohm.   f .   +-30  KHz 
tap  input 


132^  variation 


4  .. 
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2  ' 

1  . 


.04 
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Tap  Output  Amplitude 
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FIGURE  3.23  Non  Uniformity  Data,   f    »2.5  MHz 

samp 

R.   »  1  Meg-ohm.   f .    .=100  KHz 
tap  input 
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Tap  #1  Output 


Tap  #5  Output 


Tap  #6  Output 


Tap  #11  Output 


FIGURE  3.26  Pictorial  Presentation  of  TAD- 12  Performance. 

100  Kilo-ohms. 

50  KHz. 


f         »  3  MHz;     R. 
samp  tap 

V.        K»   .6  volt  peak  to  peak;      f . 


input  r  r  input 

Scale:     H  ■  10  microsec./div,   V  -   .2  volt/div 
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Tap  #1  Output 


Tap  #4  Output 


Tap  #2  Output 


Tap  #3  Output 


FIGURE  3.27  Pictorial  Presentation  of  Ret icon  TAD- 12 

Performance,  f    »  3-5  MHz.  R   -  100  Kilo-ohms 
samp  wap 

7input"  *6  V0lt  Peak  t0  Pe8jC*   f input"  3°  KHZ* 
Scale:  H  -  10  microsec./div;  V   -  .2  Volt/div 
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Tap  #1  Output 


Tap  #2  Output 


Tap  #3  Output 


Tap  #12  Output 


FIGURE  3.28  Pictorial  Presentation  of  Reticon  TAD-12. 

Performance,  f    -  4  MHz.  R.   =100  Kohms. 
samp  tap 

V    .-  .6  V  peak  to  peak;   f .    .-  50  KHz. 
Scale:  H  »  10  micro sec./div;   7  =  .2  V/div. 
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Tap  #1  Output 


Tap  #9  Output 


FIGURE  3.29  Reticon  TAD-12  Performance. 

f    ■  4.5  MHz.  R.   -  100  Kilo-ohms, 
samp  tap 

V.    .-  .6  V  peak  to  peak.  f,    .-  50  KHz, 
input      *      r  input 

Scale:  H  =  10  microsec./div;   V  *  .2  V/div. 


Tap  #1  Output 


Tap  #2  Output 


FIGURE  3.30  Ret icon  TAD-12  Performance. 

f    =5  MHz.  R.   =100  Kilo-ohms, 
samp  tap 

V.    .«  .6  V  peak  to  peak,  f .    ,-  50  KHz. 
input       c  *  input 

Scale:  H  -  10  microsec./div;   V  ■  .2  V/div, 


70 


Input 


Tap  #1   Output 


FIGURE  3.31     Reticon  TAD-12  Performance. 

f         -  2.5  MHz.     R.      =»  10  Kilo-ohms, 
samp  tap 

V.        .-   .6  V  peak  to  peak.     f.        .-  100  KHz. 
input  v  r  input 

H  ■  5  micro  sec. /div;      7  =   .2  V/div 


Upper  Scale: 

Lower  Scale:  H  =  5  microsec./div; 


V  -  .1  V/div 


Input 


Tap  #1  Output 


FIGURE  3.32  Reticon  TAD-12  Performance 


samp 


R.   =10  Kilo-ohms, 
tap 


input 
Upper  Scale:  H 


2.5  MHz. 
.6  V  peak  to  peak,  f 


.=  1  MHz. 
input 

.5  microsec./div;   V  =  .2  V/div 


Lower  Scale:   H  -   .5  microsec./div;      V  =   .1    V/div 
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Input 


Tap  #6  Output 


FIGURE  3.33  Reticon  TAD-12  Performance. 


f    -  2.5  MHz.  R.   =  10  Kilo-ohms, 
samp  tap 

V.    .»  .6  V  peak  to  peak. 
Upper  Scale:  H 
Lower  Scale :  H 


f .    .»  300  KHz. 
input 

1  microsec./div;   V  *  .2  V/div 
1  microsec./div;  V  ■  .1  V/div 


Input 


Tap  #6  Output 


FIGURE  3.34  Reticon  TAD-12  Performance. 


f         -  2.5  MHz; 
samp 

V.        A=*   .6  V  peak  to  peak,      f 


R.      =10  Kilo-ohms, 
tap 


1  MHz 


input       *       *  input 

Upper  Scale:  H  =  .5  microsec./div;  V  »  .2  V/div 

Lower  Seale:  H  »  .5  microsec./div;   V  =  .1  V/div 
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Input 


Tap  #12  Output 


FIGURE  3.35  Ret icon  TAD- 12  Performance. 


f    =  2.5  MHz. 
samp 

V. x"  .6  V  peak  to  peak,  f 


R.   =10  Kilo-ohms, 
tap 


input 
Upper  Scale:  H 

Lower  Scale:  H  ■  5  microsec./div; 


»  100  KHz 
input 

5  microsec./div;   V  =  .2  V/div 


V  *  .1  V/div 


Input 


Tap  #12  Output 


FIGURE  3.36  Ret icon  TAD- 12  Performance. 

fsamp-  2'5  mz-     Rtap-  10  Kilo-ohms. 
Vut3  .6  V  peak  to  peak,  f.^-  300  KHz 
Upper  Scale:  H  -  2  microsec./div;   V  =  .2  V/div 
Lower  Scale:  H  =  2  microsec./div;   V  ■  .1  v/div 
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Input 


Tap  #3  Output 


FIGURE  3.37  Ret icon  TAD- 12  Performance. 


f    -3  MHz.  R.   =10  Kilo-ohms, 
samp  tap 

V.    .-  .6  V  peak  to  peak;   f .    .-  200  KHz 
input       r  *  input 

Upper  Scale:  H  =  2  microsec./div;   V  -   .5  V/div 
Lower  Scale:  H  ■  2  microsec./div;   V  *  .2  v/div 


Input 


Tap  #2  Output 


FIGURE  3.38  Ret icon  TAD- 12  Performance, 
f 


samp 

r. 
mpu 

Scale: 


input 


3.5  MHz;  R.   =  10  Kilo-ohms, 
tap 

.6  V   peak  to  peak;   f .    -  200  KHz 
H  ■  2  microsec./div;   V  =  .5  V/div 
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Input 


Tap  #2  Output 


FIGURE  3.39  Reticon  TAD- 12  Performance. 


samp 
input 


R.    =10  Kilo-ohms 
tap 


5  MHz. 
.6  7  peak  to  peak. 
Upper  Scale:  H  =  5  microsec./div; 
Lower  Scale:  H  »  5  microsec./div; 


input 


200  KHz 
7  -  .2  7/div 
7  -  .01  7/div 


Input 


Tap  #1  Output 


FIGURE  3.40  Reticon  TAD- 12  Performance. 

f    =5  MHz.  R.   -  10  Kilo-ohms 
samp  tap 

7     =.6  7  peak  to  peak; 

H  =  2  microsec./div; 


f .        ,=  200  KHz 
input 


Upper  Scale: 
Lower  Scale: 


H  -  2  microsec./div; 


7 
7 


.2  7/div 
.1    7/div 
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C.   TIME  DOMAIN  EVALUATION 

The  measurement  procedure  is  to  use  a  periodic  input 
pulse.   Its  period  is  at  least  twice  the  total  delay.   The 
width  of  the  pulse  must  satisfy  the  condition: 

2   of  clock  period  _<  width  of  pulse  <  1  clock  period. 

The  pulse  amplitude  should  be  within  the  dynamic  range.   At 
the  output  of  the  summer  operational  amplifier,  one  can  ob- 
serve twelve  output  amplitudes.   An  analog  sample  and  hold 
circuit  is  used  to  measure  each  output  amplitude. 

Figures  3-45  through  3-47  illustrate  the  impulse  response 
of  TAD-12  performance  at  different  sampling  frequencies: 

fs  =  50  KHz,  1.25  MHz,  2.5  MHz 

and  tap  resistances: 

Rtap  =  10  Kfi'  10°  Kfi>  X  m 
connected  to  tap  #1,  5,  6. 
1 .   Non- Uniformity 

The  non-uniformity  data  of  tap  outputs  for  the  impulse 
input  is  shown  in  Figures  3-48,  3-49,  using  10  K^  tap  resist- 
ances at  different  sampling  frequencies: 

fs  =  25  KHz,  50  KHz,  500  KHz,  1.25  MHz. 
The  variation  of  the  tap  output  amplitudes  is  the 
following : 

-  4.21,  at  fs  =  25  KHz 

-  1.8%,  at  £   =  50  KHz 

-  4.1%,  at  f  =  500  KHZ 

-  10.471,  at  fs  =  1.25  MHz. 
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Figures  3-50,  3-51  illustrate  the  output  of  12  taps 
using  R    =  10  Kft  at  sampling  frequencies: 

fs  =  50  KHz,  1.25  MHz. 

At  lower  sampling  frequency,  up  to  f  =  1.25  MHz,  the 
non-uniformity  of  tap  output  is  less  than  10%;  when  the  samp- 
ling is  increased,  the  uniformity  is  deteriorated  seriously. 
At  f  =  2.5  MHz,  the  non-uniformity  of  different  taps  is 
higher  than  20%. 

The  interaction  between  adjacent  taps  can  also  affect 
the  tap  uniformity  at  high  sampling  frequencies.  Figure 
3-45  (c) ,  (d) ,  (e) ,  (f)  shows  these  effects.   For  example, 
in  Figure  3-45  (d)  at  f  =  1.25  MHz,  the  tap  #1  output  can 
have  an  effect  on  tap  #2  before  its  amplitude  reaches  the 
zero  value.   Another  example,  in  Figure  3-45  (f) ,  at  f  of 
2.5  MHz,  the  effect  of  tap  #1  output  can  extend  to  tap  #2 
and  tap  #3  before  its  amplitude  goes  below  the  fixed  pattern 
noise. 

2 .   Temporal  Noise 

The  noise  of  tap  outputs  can  be  grouped  into  two 
categories : 

-  fixed  pattern  noise 

-  random  noise. 

The  fixed  pattern  noise  is  periodic  with  a  period 
related  to  the  sampling  frequency  used.   At  f   of  50  KHz, 
this  period  is  about  6.67  microseconds.   At  f  of  500  KHz, 
it  is  about  80  microseconds,  but  the  number  of  periods  is 
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constant  between  tap  outputs.   For  example,  there  are  three 
periods  of  fixed  pattern  noise  between  tap  #1  output  and 
tap  #5  output.   This  is  illustrated  in  Figures  3-45  (a), 
(c) ,  3-46  (a),  (c) ,  3-47  (a),  (c) .   The  peak  to  peak  ampli- 
tude of  this  noise  is  about  40  mV  for  the  10  Kft  load  resistors 
and  about  25  mV  for  the  100  Kft  load. 

The  random  noise  varies  from  one  tap  to  the  other. 
It  depends  also  on  the  clock  frequency  and  tap  resistors. 
The  random  noise  affects  seriously  the  tap  output,  except 
tap  #1,  at  high  tap  resistor  value  and  high  clock  frequency. 
For  example,  using  10  KQ,   resistors,  the  amplitude  of  random 
noise  is  40  mV  at  clock  frequency  of  f  =  2.5  MHz  and  is 
100  mV  at  f  =  5  MHz.   This  is  shown  in  Figure  3-45  (d) ,  (f) . 
The  tap  #5  output  is  deteriorated  badly  at  f   =  5  MHz  for 

R^        =    100   Kft   and  R^        =   1  Mil   in   Figures    3-46    (d)    and 

tap  tap  &  *■  J 

3-47  (d) .   However,  the  noise  on  tap  #1  is  very  small  at  any 

condition.   This  is  illustrated  in  Figures  3-45  (d) ,  (f) , 

3-46  (d),  (f). 

The  jitter  amplitude  observed  at  the  top  of  tap  #1 

output  in  Figure  3-52  equals  to  80  microvolts  for  a  signal 

frequency  of  2.83  KHz,  f  =  100  KHz  and  for  R.    =  10  Kfi. 

c  cap 
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Figure  3-5  shows  the  magnitude  frequency  response 
of  Tap  #1,  Tap  #6  and  Tap  #8  outputs  using  100  Kilo-ohms 
resistors  and  1  Meg-ohm  resistors.   The  roll-off  starts  at 
about  20  KHz  input  frequency  for  the  R    =1  Meg- ohm  curve 
and  at  about  60  KHz  input  frequency  for  the  R    =  100  Kilo- 
ohms  curves. 

When  the  sampling  frequency  is  increased,  the  fre- 
quency dependence  curves  are  not  improved  in  comparison  with 
the  theoretical  curve.   Figures  3-6  and  3-7  illustrate  the 
magnitude  frequency  response  of  different  tap  outputs  at 
sampling  frequency  equal  to  2.5  MHz  using  tap  resistors: 

R^   =  10  Kilo-ohms,  100  Kilo-ohms,  1  Meg-ohm. 

tap  ° 

With  R    =  10  Kilo-ohms,  the  frequency  roll-off  is 
about  150  KHz.   When  the  signal  frequency,  f,  equals  to  200 
KHz,  the  measured  curve  is  1  db  below  the  theoretical  curve. 
At  f  =  300  KHz,  it  is  3  db  below  the  theoretical  curve.   For 
the  R    =  1  Meg-ohm  curves,  the  frequency  roll-off  begins 
about  30  KHz. 

With  Rtan  =  100  Kilo-ohms,  the  response  curve  begins 
to  roll  off  at  the  input  frequency  f  =  80  KHz.   At  f  =  100 
KHz,  the  response  curve  is  2  db  below  the  theoretical  curve. 

Figures  3-8  and  3-9  illustrate  respectively  the 
magnitude  frequency  response  of  TAD-12  output  using  10  Kilo- 
ohms  tap  resistors  at  sampling  frequency  fs  equal  to  2.75  MHz 
and  using  1  Meg-ohm  tap  resistor  at  fs  equal  to  1.25  MHz. 
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4.2$  variation 
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Tap  output  amplitude 
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1 .8$  variation 
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Volts 
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Tap  output  amplitude 
PIGTJRE  3.48  Non  Uniformity  Data.  Input  Pulse:  Amplitude^  Volt, 
Width=.2  microsec.  Tapping  Resistor=10  Kilo-ohms 

(a)  P    =  25  KHz 
s    samp 

(b)  P    =  50  KHz 
v  /  samp 


Volts 


Uumber 
of  tapsg 


(a) 


2  ' 
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4.1%  variation 


4-*-i 
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Tap  output  amplitude 
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10.47%  variation 
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FIGURE  3.49  Non  Uniformity  Data.  Tapping  Resistors-10  Kilo-ohms 

(a)  Input  Pulse:  Amplitude-"! .  1  Volt,  Width-. 6  microsec. 

P    «500  KHz 
samp 

(b)  Input  Pulse:  Amplitudes  .1  Volt,  Width=.3  microsec. 

P    =1.25  MHz 
samp 
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Output 


Input 


FIGURE  3.50  Impulse  Response  of  Reticon  TAD-12  showing 
12  Tap  Outputs.  R.   =10  Kilo-ohms. 
Sampling  Frequency  ■  50  KHz. 

Upper  Scale:  H  ■  50  microsec./div;  V»  .5  V/div 
Lower  Scale:  H  ■  50  microsec./div;  V  ■  1  V/div 


Output 


Input 


FIGURE  3.51   Impulse  Response  of  Reticon  TAD-12  Showing 


12  Tap  Outputs.  R 


tap 


10  Kilo-ohms. 


Sampling  Frequency  =1.25  MHz. 
Upper  Scale:  H  -  2  microsec./div;  V 
Lower  Scale:  H  =  2  micro sec. div/;   V 


.2  V/div 
1  V/div 
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Output 


Tap  #1  Jittering 
Waveform 


FIGURE  3.52  Upper  Trace:  Impulse  Response  of  12  Tap  Outputs 

Scale:  H  *  .1  msec./div;  V  -  .05  v/div 

Lower  Trace:  Jittering  Waveform  of  Tap  #1 

Scale:  H  =  .1  msec./div;  V  =  .02  V/div 

Sampling  Frequency  »  50  KHz 

Tapping  Resistance  »  10  Kilo-ohms. 
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IV.   INVESTIGATION  OF  FILTER  I  -  PREWHITENING  FILTER 

A.   THEORY 

Based  on  the  design  algorithm  described  in  II.D.l,  a 
computer  program  shown  in  Appendix  A  is  written  to  give  the 
impulse  response  of  the  prewhitening  filter  for  the  required 
frequency  characteristics.   Using  this  program,  a  prewhitening 
filter  is  designed  for  the  following  conditions: 


f  =  2.1  KHz 
o 


f„  =  1/2  f 
p        s 


=  625  KHz 


R  =  0.1  Slope  A  =  1 

Length  of  filter:   M  =  11  taps 


H(f» 


R 

I  | 

i  I  

f.  f      f  normalized 

0  p 

The  tapping  coefficients  improved  by  Hamming  window  are 
listed  in  Table  IV. I. 

The  magnitude  frequency  response  of  the  theoretical  pre^ 
whitening  filter  is  shown  in  Figure  4-1.   It  is  plotted  in 
linear  scale  on  Y-axis  versus  normalized  frequencies, 
f/f  on  X-axis. 
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TABLE  IV. I   Theoretical  Tapping  Coefficients  §  Resistors 


k 

Ak 

R,  (Kilo-ohms) 

A,  Measured 
k 

-5 

-0.00091 

-16478.516 

0.02564 

-4 

0.0 

open 

0.15385 

-3 

-0.01219 

-1225.517 

0.05128 

-2 

0.0 

open 

0.0 

-1 

-0.21407 

-65.005 

-0.23076 

0 

1. 

10. 

1. 

1 

-0.21427 

-65.005 

-0.19231 

2 

0.0 

open 

0.02564 

3 

-0.01219 

-1225.517 

0.02564 

4 

0.0 

open 

0.02566 

5 

-0.00091 

-16478.516 

-0.07692 

B.   EXPERIMENT 

1.   Tapping  Resistors  According  to  Theoretical 
Design  with  Rmin  -  10  Kilo-ohms 

A  minimum  value  of  10  Kilo-ohms  is  set  for  the  tapping 

resistances.   Appendix  F  indicates  the  following  expression 

for  Rk, 

Rv  =  (lOK  +  r)  C^2*D  -  r 


where  r  is  the  internal  output  resistance.   For  the  TAD-12, 
r  is  about  5  Kilo-ohms.   The  tapping  resistors  corresponding 
to  the  designed  prewhitening  filter  are  given  in  Table  IV. 1 
and  used  to  implement  the  filter  using  the  11  taps  of 
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Reticon  TAD-12.   Since  their  theoretical  coefficients  A, 

k 

for  k  higher  than  5  have  negligible  values,  the  implementa- 
tion of  the  filter  on  a  longer  delay  line,  number  of  taps 
higher  than  11,  is  not  necessary. 

The  measured  impulse  response  of  this  filter  is 
illustrated  in  Figure  4-2  at  three  different  sampling  fre- 
quencies : 

fs  =  50  KHz,  500  KHz,  1.25  MHz 

The  pulse  inputs   used  have  the  amplitude  and  width  listed 
in  Table  IV. II. 

TABLE  IV. II   Pulse  Inputs 


fs 

Amplitude 

Width 

50  KHz 

1.5  V 

8  microsec. 

500  KHz 

1.5  V 

.  8  microsec. 

1.25  MHz 

1.5  V 

.  3  microsec. 

Using  the  sample  and  hold  circuit,  the  impulse 
response,  A,  is  measured  and  the  result  is  listed  in  Table 
IV. 1.   The  A,  measured  and  the  theoretical  A,  are  quite 
different  for  some  taps.   The  maximum  deviation  occurs  be- 
tween tap  #7  and  #8  outputs.   Its  peak  amplitude  is  about 
30%  of  the  output  at  tap  #6. 

The  theoretical  magnitude  frequency  response  and 
measured  response  curves  for  three  different  sampling  fre- 
quencies are  plotted  on  Figure  4-3.   The  response  curves  for 
f   of  1.25  MHz  and  500  KHz  show  some  wiggling  shape  due 
probably  to  the  non-uniformity  of  tap  output  and  the  clock 
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noise  which  cause  the  unsymmetrical  tapping  coefficients. 
The  response  curve  for  f  equal  to  50  KHz  is  rather  satis- 
factory if  it  is  multiplied  by  a  factor  of  1.07. 

2 .   Tapping  Resistors  Adjusted  to  Yield 

Better  Impulse  Response.  Rmin  =  10  Kilo-ohms 

Since  only  three  taps  have  strong  contribution  to 

the  filter,  the  adjustment  is  done  only  on  the  taps  #5,  #6, 

#7  by  using  the  sample  and  hold  circuit.   The  other  taps  are 

left  open.   The  tapping  resistors  used  are  listed  in  Table 

IV. III. 

TABLE  IV. III.   Tapping  Resistors 

k        Av  Rk  (Kilo-ohms) 

K      fs  =  50  KHz   fs  =  500  KHz   fs  =  1.25  MHz 

-1     -0.21427       28  50  52 

0  1.  10  10  10 

1  -0.21427       28  35  100 

The  amplitude  and  width  of  inputs  pulse  are  listed  in 
Table  IV. II. 

The  impulse  response  of  this  filter  is  illustrated  in 
Figure  4-4,  at  three  different  sampling  frequencies. 

The  theoretical  and  measured  curves  of  prewhitening 
filter  magnitude  frequency  response  are  shown  in  Figure  4-5. 
The  response  curves  for  fg  equal  to  50  KHz,  500  KHz  and  1.25  MHz 
still  show  some  wiggling.   The  gain  of  these  curves  are  lower 
than  the  last  curves  in  Figure  4-3.   If  these  three  curves 
are  multiplied  by  a  factor  of  1.22,  they  are  somewhat  satis- 
factory in  comparison  with  the  theoretical  curve. 
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3.   Tapping  Resistors  Adjusted  to  Yield 

Better  Impulse  Response.   Rmin  =  50  ohms 

A  minimum  value  of  50  ohms  is  set  for  the  tap  #6. 

The  taps  #5  and  #7  are  adjusted  by  using  sample  and  hold 

circuit,  to  yield  better  impulse  response.   The  other  taps 

are  left  open.   The  tapping  resistors  used  are  listed  in 

Table  IV. IV. 

TABLE  IV. IV.   Tapping  Resistors 

k      A]<  Rfc  (Kilo-ohms) 

fs  =  50  KHz    fs  =  500  KHz    fs  =  1.25  MHz 


-1 

-0.21427 

4.8 

2.8 

3.8 

0 

1. 

0.05 

0.05 

0.05 

1 

-0.21427 

2.8 

1.8 

11.5 

The  pulse  inputs  used  have  the  amplitudes  and  widths 
listed  in  Table  IV. II. 

Figure  4-6  illustrates  the  impulse  response  of  the 
filter  at  sampling  frequencies: 

fs  =  50  KHz,  500  KHz,  1.25  MHz 

The  theoretical  and  measured  curves  of  prewhitening 
filter  magnitude  frequency  response  are  shown  in  Figure  4-7. 
The  wiggling  of  the  three  measured  curves  is  much  reduced  in 
comparison  with  the  previous  curves.   If  the  experimental 
curve  f  =  1-25  MHz  is  multiplied  by  a  factor  of  1.2  and  curves 
f   =  50  KHz,  500  KHz  multiplied  by  a  factor  of  1.35,  the 
agreement  is  somewhat  satisfactory. 
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4.   Tapping  Resistors  Obtained  by  Experimental 

Cut  $  Try  Based  on  Frequency  Response  Adjustment 

The  procedure  consists  of  the  use  of  a  chirped  wave- 
form from  DC  to  1  MHz  as  the  input  of  the  filter,  and  adjust 
the  displayed  frequency  response  to  get  the  minimum  wiggling 
by  adjusting  the  tapping  resistances.   The  result  obtained 
is  shown  in  Figure  4-9.   The  wiggling  of  the  three  response 
curves  is  much  reduced.   Raising  these  curves  by  a  factor  of 
1.27  can  improve  the  agreement  between  the  theoretical  and 
measured  curves. 

The  impulse  response  corresponding  to  these  curves 
is  unsymmetrical  and  has  the  values  listed  in  Table  IV. V. 

TABLE  IV.V   Tapping  Coefficients 

k  A^       Rk  (Kilo-ohms) 


-1        -0.1694       24.812 

0  1.  .05 

1  -1  .05 

The  amplitude  and  width  of  pulse  inputs  used  are 
listed  in  Table  IV. II. 
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FIGURE  4.1  Theoretical  Amplitude  Frequency  Response  of  Prewhitening 
Filter.   Length  -11.  Hamming  window. 
X  -  Scale  =  .1  units/inch  for  normalized  frequency 
Y.  -  Scale  =»  .2  units/inch  for  Amplitude  frequency  response, 
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V.   INVESTIGATION  OF  FILTER  II  -  DEWHITENING  FILTER 

A.   THEORY 

The  design  algorithm  of  the  dewhitening  filter  is 
described  in  II. D. 2.   The  computer  program  in  Appendix  B  is 
used  to  calculate  the  impulse  response  of  the  filter  from 
the  given  frequency  response.   A  dewhitening  filter  is 
designed,  based  on  the  impulse  response  of  the  prewhitening 
filter  of  paragraph  IV. A.   The  frequency  characteristics 
are  the  following: 


£n  =  2.1  KHz 


fp  =  |  fs  =  625  KHz 


Length  of  filter  M  =  11  taps 


K-  10 


l*«u 


H 


f     f  normalized 
P 


The  tapping  coefficients  are  improved  by  a  Hamming  window 
and  by  a  factor  of  1.13  and  are  listed  in  Table  V.I. 

The  magnitude  frequency  response  of  the  theoretical  pre- 
whitening filter  is  shown  in  Figure  5.1.   It  is  plotted  in 
linear  scale  on  Y-axis  versus  normalized  frequencies, 
f/fs  on  X-axis. 
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TABLE  V.I.  Theoretical  Tapping  Coefficients  §  Resistors 

k    Afc  theor.      R^  (Kilo-ohms)      R^  (Kilo-ohms) 

with  Rmin  =  K  ohms   with  Rmin  =  50  ohms 


-5  0.00113  13269.336  4464.026 

-4  0.00632  2368.417  794.050 

-3  0.01156  1292.577  431.851 

-2  0.04846  304.661  99.252 

-1  0.19681  71.215  20.659 

0  1.  10.  .05 

1  0.19681  71.215  20.659 

2  0.04844  304.661  99.252 

3  0.01156  1292.577  431.851 

4  0.00632  2368.417         .   794.050 

5  0.00113  13269.336  4464.026 

B.   EXPERIMENT 

1.   Tapping  Resistors  According  to  Theoretical 
Design  with  Rmin  =  10  Kilo-ohms 

A  minimum  value  of  10  Kilo-ohms  is  set  for  the 
tapping  resistors.   Using  the  expression  of  R^  given  in 
Appendix  F,  the  tapping  resistances  are  calculated  and 
listed  in  Table  V.I  in  the  third  column.   They  are  used  to 
implement  the  dewhitening  filter  using  11  taps  of  Reticon 
TAD- 12. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  5-2  at  three  different  sampling  frequencies: 

fs  =  50  KHz,  500  KHz,  1.25  MHz. 
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5 

0.00113 

4 

0.00632 

3 

0.01156 

2 

0.04844 

1 

0.19681 

0 

1. 

1 

0.19681 

2 

0.04844 

3 

0.01156 

4 

0.00632 

5 

0.00113 

TABLE  V.II   Theoretical  §  Measured  Coefficients 


k   A^  theor.  A^  Measured 

fs  =  50  KHz    fs  =  500  KHz    fs  =  1.25  MHz 


0.16216  0.10869  0.11667 

0.16216  0.17391  0.20830 

0.16216  0.15217  0.25 

0.2027  0.19565  0.2708 

0.39119  0.41304  0.41667 

1.  1.  1. 

0.39189  0.39130  0.41667 

0.18918  0.13043  0.0 

0.16216  0.13043  0.16667 

0.17567  0.13043  0.16667 

0.16216  0.10869  0.16667 


The  impulse  inputs  are  noted  in  Table  IV. II  for 
different  sampling  frequencies.   The  tapping  coefficients 
are  measured  with  the  sample  and  hold  circuit  and  the  results 
are  listed  in  Table  V.II.   The  theoretical  and  measured  A^ 
are  quite  different  for  some  taps  due  to  the  fixed  pattern 
noise  and  non-uniformity  of  tap  outputs.   The  maximum  noise 
amplitude  is  observed  before  the  tap  #1  output  and  has 
amplitude  of  about  30%  of  A0. 

The  theoretical  magnitude  frequency  response  and 
measured  response  for  three  different  sampling  frequencies 
are  plotted  in  Figure  5-3.   The  response  curves  for  fs 
of  500  KHz  and  1.25  MHz  are  somewhat  in  agreement  with  the 
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theoretical  curve.   Their  difference  is  about  12%  for  f 

less  than  0.2  f  and  it  becomes  worse  than  40%  for  input 

signal  frequency  higher  than  0.25  f  .   For  f  of  50  KHz, 

the  result  is  quite  satisfactory. 

2.   Tapping  Resistors  Adjusted  to  Yield  Better 
Impulse  Response.   Rm-jn  =  10  Kilo-ohms 

Since  the  coefficients,  A,  for  k  greater  than  3, 
have  negligible  values  in  comparison  with  the  noise,  the 
adjustment  is  done  only  on  seven  tap  outputs.   The  adjust- 
ment is  quite  delicate  for  the  tap  outputs  which  are  greater 
than  or  equal  to  5,  because  of  the  interaction  between  them. 
Experimental  tap  resistances  are  listed  in  Table  V.III.   A 
slight  modification  has  to  be  done  on  the  values  of  some 
tap  resistors  in  order  to  get  better  impulse  responses  for 
different  sampling  frequencies.   The  sample  and  hold  circuit 
is  used  to  make  the  adjustment. 

The  amplitude  and  width  of  inputs  pulse  used  are 
listed  in  Table  IV. II. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  5-4  at  three  different  sampling  frequencies.   The 
maximum  noise  occurs  near  tap  #1  output;  its  peak  amplitude 
is  about  19%  of  AQ. 

The  theoretical  and  measured  curves  of  dewhitening 

filter  magnitude  frequency  response  are  shown  in  Figure  5-5. 

The  wiggling  of  the  response  curves  is  reduced  in  comparison 

with  the  curves  in  Figure  5-3,  but  the  difference  between  the 

theoretical  and  measured  curves  (f  =1.25  MHz  and  500  KHz) 

v  s 
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is  still  high,  greater  than  301,  for  the  input  signal 
frequency  f  higher  than  0.25  f  .   At  low  sampling  frequency, 
50  KHz,  the  result  is  quite  satisfactory. 


TABLE  V.III.   Adjusted  Coefficients  §  Resistors 


Ak  Rk 

f  =50  KHz     f   =500  KHz  f   =  1.25MHz      f   =50  KHz     f   =500  KHz  f   =1.25MHz 
s  s  s  s  s  s 


5 

0.09836 

0.02083 

0.125 

4 

0.10655 

0.02083 

0.1 

3 

0.10655 

0.04167 

0.075 

2 

0.09836 

0.0625 

0.075 

1 

0.18852 

0.20833 

0.175 

0 

1. 

1. 

1. 

1 

0.18852 

0.20833 

0.2 

2 

0.09836 

0.0625 

0.1 

3 

0.10655 

0.02083 

0.075 

4 

0.11675 

0.02083 

0.075 

5 

0.10655 

0.02083 

0.075 

13269. 
2368. 
1292.5 
300. 
110. 
5. 
110 
311. 
1292.5 
2368. 
13269. 


13269. 

2368. 

1292.5 

300. 

109. 

9.5 

110. 

311. 

1295.5 

2368. 

13269. 


13269. 

2368. 

1292.5 

300. 

110. 

5. 

110. 

320. 

1292.5 

2368. 

13269. 


3 .   Tapping  Resistors  According  to  Theoretical  Design. 
Rniin  =  50  ohms. 

A  minimum  value  of  50  ohms  is  set  for  the  tap  #6.   Using 
the  formula  described  in  Appendix  F,  one  can  calculate  the 
theoretical  tapping  resistances  whose  values  are  listed  in  the 
fourth  column  of  Table  V.I. 

The  input  pulses  have  the  amplitude  and  width  listed 
in  Table  IV. II.   Figure  5-6  illustrates  the  impulse  response 
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of  the  filter  at  sampling  frequencies: 

f   =  50  KHz,  500  KHz,  1.25  MHz 

s  ' 

The  measured  coefficients  A,  are  listed  in  Table  V. IV 


TABLE  V.IV.   Measured  Coefficients  (R  .  =50  ohms) 

v  mm 


k 

f  =50  KHz 
s 

r  Ak 

fs=500  KHz 

f  =1.2  5  MHz 
s 

-5 

0.11111 

0.06667 

0.22826 

-4 

0.11111 

0.06667 

0.08695 

-3 

0.11111 

0.06667 

0.09782 

-2 

0.16667 

0.13333 

0.15217 

-1 

0.27778 

0.26667 

0.34782 

0 

1. 

1. 

1. 

1 

0.22223 

0.26667 

0.48913 

2 

0.13888 

0.13333 

0.20652 

3 

0.11111 

0.06667 

0.13043 

4 

0.11111 

0.06667 

0.04347 

5 

0.11111 

0.06667 

0.07608 

The  maximum  peak  amplitude  of  noise  measured  is  about 

131  of  A  for  f  equal  to  1.25  MHz  and  500  KHz  and  191  of  A 
o      s   n  o 

for  f   equal  to  100  KHz. 

The  theoretical  and  measured  curves  of  dewhitening 
filter  magnitude  frequency  response  are  shown  in  Figure  5-7. 
The  response  curve  for  f   equal  to  50  KHz  is  similar  to  the 
theoretical  one.   The  curves  f   equal  to  1.25  MHz  and  500  KHz 
have  some  wiggling  effect  for  input  frequencies  higher  than 

.29  fs.  112 


Required  DWF 


X  X  Designed  DWF  improved 
by  Hamming  window  and 
a  factor  of  1.1 3 


fe> 


FIGURE  5.1  Amplitude  Frequency  Response  of  Dewhitening  Filter, 
Length  -  11  /  Hamming  window 

X  -  Scale  -  .1  units/inch  for  normalized  frequency 
Y  -  Scale  ■  .5  units/inch  for  amplitude  frequency 
response. 
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VI.   INVESTIGATION  OF  FILTER  III  -  BANDPASS  FILTER 

A.   THEORY 

Based  upon  the  design  algorithm  described  in  II. D. 3,  a 
computer  program  shown  in  Appendix  C  is  written  to  calculate 
the  impulse  response  of  the  bandpass  filter  based  on  the 
required  frequency  response.   Using  this  program,  a  bandpass 
filter  is  designed  for  the  following  conditions: 


f2  =  .1996 


Center  frequency  f 

Amplitude  A 

Number  of  taps  M 

|H(f)  | 


.2 
1. 
11 


f2  =  .2004 


Q  =  250 


A  •• 


i »■ 


h  -2  h 


f  normalized 


The  tapping  coefficients  improved  by  Hamming  window  are 
listed  in  Table  VI. I. 

The  magnitude  frequency  response  of  the  theoretical  band- 
pass filter  is  shown  in  Figure  6-1,  in  linear  scale,  and  in 
Figure  6-2,  in  logarithmic  scale  versus  normalized  frequencies, 
f/f  . 
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TABLE  VI. I.   Theoretical  Coefficients  and  Resistors 


k 

Ak 

R,  (R    .    =10  k) 
kK  mm 

k v  mm 

-5 

0.09863 

147.083 

46.201 

-4 

0.07378 

198.307 

63.447 

-3 

-0.38389 

34.074 

8.155 

-2 

-0.59148 

20.360 

3.538 

-1 

0.28647 

47.362 

12.628 

0 

1. 

10. 

.05 

1 

0.28647 

47.362 

12.628 

2 

-0.59148 

20.360 

3.538 

3 

-0.38389 

34.074 

8.155 

4 

0.07378 

198.307 

63.447 

5 

0.09863 

147.083 

46.201 

B.   EXPERIMENT 

1.   Tapping  Resistors  According  to  Theoretical 
Design  with  Rm-in  =  10  K  ohms. 

A  minimum  value  of  10  Kilo-ohms  is  set  for  the 
tapping  resistors.   The  expression  of  R,  given  in  Appendix  F 
is  used  to  calculate  the  tapping  resistances.   These  values 
are  listed  in  Table  VI . I  in  the  third  column  and  are  used  to 
implement  the  bandpass  filter  using  the  11  taps  of  Reticon 
TAD- 12. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  6-3  at  three  different  sampling  frequencies: 

f   =  50  KHz,  500  KHz,  1.25  MHz 

s         ' 
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5 

0.09863 

4 

0.07378 

3 

-0.38389 

2 

-0.59148 

1 

0.28647 

0 

1. 

1 

0.28647 

2 

-0.59148 

3 

-0.38389 

4 

0.07378 

5 

0.09863 

TABLE  VI. II.   Theoretical  $  Measured  Coefficients 

k   A,  Theor.  A^  Measured 

f  =50  KHz      f  =500  KHz     f  =1.25  MHz 

55 S S 

0.34444  0.30508  0.34483 

0.33333  0.209339  0.06896 

-0.5  -0.66102  -0.86306 

-0.6  -0.66102  -0.51724 

0.48888  0.49152  0.72413 

1.  1.  0.86206 

0.44444  0.54237  0.4 

-0.65555  -0.76271  -1. 

-0.48888  -0.71186  -0.34482 

0.31111  0.13559  0.51724 

0.26666  0.13559  0.24138 

The  input  pulses  used  are  listed  in  Table  IV. II  for 
different  sampling  frequencies.   The  tapping  coefficients 
measured  with  the  sample  and  hold  circuit  are  shown  in  Table 
VI. II.  The  theoretical  and  measured  A,  are  quite  different 
for  some  taps  due  to  the  fixed  pattern  noise  and  the  non- 
uniformity  of  the  tap  outputs.   The  maximum  noise  observed 
occurs  between  tap  #3  and  #4  outputs.   Its  peak  amplitude  is 
about  30%  of  AQ. 

The  theoretical  magnitude  frequency  response  and 
measured  response  curves  for  three  different  sampling  fre- 
quencies are  plotted  in  Figure  6-4.   The  main  lobes  of  measured 
response  curves  are  rather  similar  to  the  theoretical  one. 
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But  there  are  some  differences  at  the  side  lobes  which  have 
their  peak  magnitudes  as  follows: 

-  15%  of  maximum  amplitude  for  f  =  50  KHz 

r  s 

-  35%  of  maximum  amplitude  for  f  =  500  KHz 

-  51%  of  maximum  amplitude  for  f  =1.25  MHz. 
There  are  also  some  wiggles  for  input  frequencies  less  than 

.15  f  - 

s 

2.   Tapping  Resistors  Adjusted  to  Yield  Better 
Impulse  Response,  Rm-jn  =  10  Kilo-ohms 

Since  for  k  greater  than  3,  the  coefficients  A,  have 
negligible  values  in  comparison  with  the  noise,  the  adjustment 
is  done  only  on  7  tap  outputs.   Obtaining  the  exact  values  of 
the  theoretical  design  is  very  delicate  because  of  the  inter- 
action between  tap  outputs.   The  adjusted  coefficients  and 
the  corresponding  tap  resistances  are  listed  in  Table  VI. III. 
Some  modification  has  to  be  done  on  the  values  of  some  tap 
resistors  in  order  to  get  better  impulse  response  for  differ- 
ent sampling  frequencies. 

The  amplitude  and  width  of  input  pulses  are  listed 
in  Table  IV. II. 

The  impulse  response  of  this  filter  is  illustrated 
in  Figure  6-5  at  three  different  typical  sampling  frequencies. 
The  maximum  noise  observed  has  its  peak  amplitude  about  19% 


of  A0. 

The  theoretical  and  measured  curves  of  bandpass  filter 
magnitude  frequency  response  are  shown  in  Figure  6-6.   The 
wiggling  of  the  frequency  response  at  input  frequencies  less 
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than  .15  f   is  much  reduced  for  f  =  50  KHz,  but  the  side 
s  s         ' 

lobe  is  still  high  at  15%  of  the  maximum  amplitude.   For  f  = 
1.25  MHz  and  f  =  500  KHz,  the  peak  amplitudes  of  side  lobes 
are  increased  to  29%  of  maximum  amplitude. 


TABLE  VI. III.   Adjusted  Coefficients  §  Resistors  (Rmin~  10  k) 

k  Afc  Rk  (Kilo-ohms) 

f  =50KHz  f  =500KHz   f  =1.25MHz  f  =50KHZ  f  =500KHz  f  =1.25MHz 
s         s  s         s        s         s 


-5 

.07246 

.01587 

.07317 

147. 

147. 

147. 

-4 

.08695 

.04762 

.04878 

198. 

198. 

198. 

-3 

-.20289 

-.55555 

-.34146 

48.07 

40. 

60.11 

-2 

-.62318 

-.62492 

-.53658 

23. 

20. 

16. 

-1 

.28985 

.30156 

.26829 

47. 

75. 

210. 

0 

1. 

1. 

1. 

8. 

8. 

6. 

1 

.28985 

.30156 

.26829 

47. 

75. 

9.8 

2 

-.62318 

-.63492 

-.53658 

23. 

24.8 

8.5 

3 

-.20289 

-.55555 

-.34146 

48.07 

40. 

23. 

4 

.08695 

.04762 

.07317 

198. 

198. 

198. 

5 

.07246 

.01587 

-.04878 

147. 

147. 

147. 

3.   Tapping  Resistors  According  to  Theoretical 
Design.   Rmin  =  50  ohms 

A  minimum  value  of  50  ohms  is  set  for  the  tapping 
resistance.   Using  the  formula  described  in  Appendix  F,  the 
tapping  resistances  are  calculated  and  have  the  values  listed 
in  Table  VI . I ,  fourth  column. 

The  pulse  inputs  used  have  the  amplitude  and  width 
listed  in  Table  IV. II.   Figure  6-7  illustrates  the  impulse 
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the  impulse  response  of  the  filter  at  sampling  frequencies 

f   =  50  KHz,  500  KHz,  1.25  MHZ 
s 

The  measured  coefficients  A,  are  listed  in  Table  VI. IV. 


TABLE  VI. IV.   Measured  Coefficients  (R  .  =  50  ohms) 


k 

fs=50  KHz 

Ak 

f  =500  KHz 
s 

f  =1.25  MHz 
s 

-5 

.03211 

.03137 

.82810 

-4 

.05504 

.07058 

.82062 

-3 

-.37155 

-.27451 

.76083 

-2 

-.58256 

-.40392 

.74439 

-1 

.24312 

.22745 

.84603 

0 

1. 

1. 

1. 

1 

.24312 

.24314 

.84603 

2 

-.58256 

-.34509 

.73243 

3 

-.37155 

-.38823 

.73094 

4 

.05504 

.04313 

.83408 

5 

.05504 

.09803 

.81464 

The  maximum  peak  amplitude  of  noise  measured  is  about  181  of  A0 

The  theoretical  and  measured  curves  of  bandpass  filter 
magnitude  frequency  response  are  shown  in  Figure  6-8.   The 
wiggling  effect  of  the  response  curves  are  strongly  diminished 
in  comparison  with  previous  curves  in  Figures  6-4  §  6-6.   The 
amplitudes  of  side  lobes  are  also  reduced  and  have  the  follow- 
ing values: 
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7.5%   of  maximum  amplitude  for  f  =  50  KHz 

-  11%  of  maximum  amplitudes  for  f  =  500  KHz 

-  12.5%   of  maximum  amplitude  for  £  =  1.25  MHz. 

4.   Tapping  Resistor  Adjusted  to  Yield  Better 
Impulse  Response.   Rmin  =  50  ohms 

a.   Driving  Circuit  Using  RCA/CMOS  Flip-Flop 

A  minimum  value  of  50  ohms  is  set  for  tap  #6. 

The  tap  outputs  are  adjusted  by  sample  and  hold  circuit  to 

yield  better  impulse  response.   The  tapping  coefficients  and 

resistances  adjusted  are  listed  in  Table  VI. V. 

TABLE  VI. V.   Adjusted  Tapping  Coefficients  §  Resistors 
Ak  Rk  (Kilo-ohms) 


f  =50KHz  f  =500KHz  f  =1.25  MHz  f  =50KHz  f  =500KHz  f  =1.25  MHz 


k 
=50KH; 


-5 

.16733 

.08532 

.01234 

46.2 

46.2 

46.2 

-4 

.06374 

.07508 

.03086 

63.45 

63.45 

63.45 

-3 

-.37450 

-.33673 

-.25308 

6.25 

4.05 

4.05 

-2 

-.60557 

-.59183 

-.59876 

1.2 

1.2 

.5 

-1 

.28287 

.23849 

.22222 

23.023 

25.023 

10. 

0 

1. 

1. 

1. 

.05 

.05 

.05 

1 

.28287 

.23849 

.23457 

6.1 

6.1 

4.1 

2 

-.60557 

-.59183 

-.59876 

1.083 

1.083 

1.08 

3 

-.37450 

-.33673 

-.23457 

4.083 

4.55 

1.083 

4 

.06374 

.07508 

.04321 

63.45 

63.45 

63.45 

5 

.16733 

.08532 

.02469 

46.2 

46.2 

46.2 

The  input  pulse  used  have  the  amplitudes  and  widths  listed  in 
Table  IV. II.   Figure  6-9  illustrates  the  impulse  response  of 
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the  filter  at  three  different  typical  sampling  frequencies. 
The  maximum  noise  has  its  peak  amplitude  about  15%  of  A  . 

The  theoretical  and  measured  curves  of  bandpass 
filter  magnitude  frequency  response  are  shown  in  Figure  6-10. 
The  wiggling  does  not  appear  at  the  main  lobe  but  the  ampli- 
tudes of  the  side  lobes  are  still  high: 

-  15%  of  maximum  amplitude  for  f  =  500  KHz 

-  10%  of  maximum  amplitude  for  f  =  500  KHz  §  1.25  MHz 
b.   Driving  Circuit  Using   INSELEX  CMOS/SOS  Flip-Flop 

A  new  circuit  with  faster  flip-flop  is  made  and  a 

minimum  value  of  50  ohms  is  set  for  tapping  resistance.   The 

adjustment  is  done  for  f  equal  to  500  KHz;  the  coefficients 

and  resistances  obtained  are  used  also  for  f  of  50  KHz  and 

s 

1.25  MHz.   Table  VI. VI  lists  these  coefficients  and  resistances. 

TABLE  VI. VI   Adjusted  Coefficients  §  Resistances 

k  A^  R^ (Kilo- ohms) 

-5  .03536  46.2 

-4  .01571  63.45 

-3  -.40078  3.8 

-2  -.61689  1.56 

-1  .29194  9.023 

0  1.  .05 

1  .29273  2.32 

2  -.61297  1.263 

3  -.40078  5.000 

4  .06287  63.45 

5  .07269  46.2 
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The  pulse  inputs  listed  in  Table  IV. II  are  used.   Figure  6-11 

illustrates  the  impulse  response  of  the  filter  at  sampling 

frequencies . 

f   =  50  KHz,  500  KHz,  1.25  MHz. 
s 

The  maximum  peak  amplitude  noise  observed  is  about  21.61  of 

v 

The  theoretical  and  measured  curves  of  bandpass 
filter  magnitude  frequency  response  are  shown  in  Figure  6-12. 
There  is  no  wiggling  at  the  main  lobe  region  and  some  improve- 
ment is  obtained  at  the  side  lobes,  which  have  the  following 
amplitudes : 

6%  of  maximum  amplitude  for  f   =  1.25  MHz 

-  10%  of  maximum  amplitude  for  f  =  500  KHz 

-  14%  of  maximum  amplitude  for  f  =  50  KHz. 

Using  the  computer  program  shown  in  Appendix  E, 
the  measured  magnitude  frequency  response  for  f     equal  to 
500  KHz  and  the  response  corresponding  to  the  measured  impulse 
response  listed  in  Table  VI. VI  are  plotted  in  Figure  6-13  with 
DB  scale  on  Y  axis  and  normalized  frequencies  on  X  axis.   The 
main  lobes  of  the  two  curves  are  the  same  but  there  are  some 
differences  at  the  side  lobe  regions.   This  shows  that  the 
technique  using  the  sample  and  hold  circuit  to  adjust  the 
tapping  resistor  to  get  better  impulse  response  is  quite 
accurate  and  very  useful  to  the  filter  implementation. 
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5.   Distortion  Measurements  of  Bandpass  Filter 

For  the  case  described  in  V.B.3  where  the  tapping 
resistances  are  set  according  to  theoretical  design  with 

R    equal  to  50  ohms,  the  second  harmonic  distortion  de- 

mm  n  ' 

fined  by: 

nR   20  1    2nd  harmonic  voltage 
°^  fundamental  voltage 

are  illustrated  respectively  in  Figures  6-14  and  6-15,  for 
sampling  frequency  of  500  KHz  and  three  different  peak  to 
peak  input  voltages: 

V.   =  3V,   0.8V,   0.1 V 
in 

From  the  Figure  6-14  where  the  input  frequency  is  equal  to 

50  KHz,  the  second  harmonic  and  clock  noise  at  85  KHz  have 

the  relative  values  expressed  in  percentage  of  fundamental 

as  shown  in  Table  VI. VI I. 

TABLE  VI. VII   Second  Harmonic  §  Clock  Noise 


Vin  ppk 

3  Volts 

.8  Volts 

.1  Volt 

2nd  Harmonic 

16.7% 

6.5% 

6.8% 

Clock  Noise 

5.8% 

14.9% 

88.8% 

When  the  amplitude  of  the  input  voltage  decreases,  the  rela- 
tive value  of  second  harmonic  decreases  but  the  relative  value 
of  clock  noise  increases.   In  the  Figure  6-14c,  one  can  observe 
many  spurious  responses  which  occur  at  63  KHz,  85  KHz,  112  KHz, 
126  KHz,  and  are  due  principally  to  the  clock  noise  and  fixed 

pattern  noise. 

I 
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FIGURE  6.1   Amplitude  Frequency  Response  of  Bandpass  Filter 
Length  =  11  /  Hamming  Window 

X  -  Scale  =  .2  units/inch  for  normalized  frequency 
Y  -  Scale  =  10  units/inch  for  amplitude  frequency 
response. 
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FIGURE  6.2  Amplitude  Frequency  Response  of  Bandpass  Filter  . 
Length  »  11  /  Hamming  Window 

X  -  Scale  *  .2  units/inch  for  normalized  frequency 
Y  -  Scale  «  10  units/inch  in  DB  for  amplitude  frequency 
response. 
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Curve  corresponding  to 
measured  impulse  response 
D  Measured  curve 


FIGURE  6.13 


Bandpass  Filter  Magnitude  Frequency  Response 
Tapping  Resistors  Adjusted  Using  Impulse  Response. 

!HnSca: 

Y  -  Scale  =  5  units/inch  in  DB 


50  ohms.   Q  =  250.  _fSam] 
le  =  .1  units/men  for  nt 


=  500  KH: 
rmalized  : 


requency 
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6.   Observations 

The  characteristic  feature  of  the  filter  realized 
by  the  electronic  tapped  analog  delay  line  is  that  one  can 
scan  the  center  frequency  of  the  bandpass  filter  to  the  left 
or  right  by  simply  decreasing  or  increasing  the  clock  fre- 
quency.  Figure  6-16  illustrates  this  feature  for  Q  =  constant 
at  two  different  sampling  frequencies: 

f   =  500  KHz,  835  KHz 


f  =500  KHz 
s 


f  =835  KHz 
s 


FIGURE  6-16. 


Bandpass  Filters  Amplitude  Frequency 
response  at  2  different  sampling 
frequencies  f  =500  KHz  d,   f  =835  KHz. 
Scale :  H  =  .  2snsec/Div  V  =s .  1  V/Div. 


Another  observation  is  that  the  transversal  filter 
implementation  is  not  very  sensitive  to  high  specifications 
for  moderate  number  of  taps.   The  tapping  coefficients  are 
almost  the  same  for  different  values  of  Q  required.   The 
impulse  response  of  bandpass  filters  for: 

fQ  =  0.2  fs     (f  =  center  frequency) 

Q  =  50,   250,   500 
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using  11  taps  is  listed  in  Table  VI. VIII,  and  using  23  taps 
in  Table  I. IX. 

TABLE  VI. VI I I   Tapping  Coefficients  (11  Taps) 


k 

Q=50 

Ak 

Q  =  250 

Q=500 

0 

1. 

1. 

1. 

1 

0.28644 

0.28647 

0.28640 

2 

-0.59140 

-0.59148 

-0.59148 

3 

-0.38382 

-0.38389 

-0.38395 

4 

0.07375 

0.7378 

0.07378 

5 

0.09857 

0.09863 

0.09863 

TABLE  VI . IX   Tapping  Coefficients  (23  Taps) 

k  Ak 
Q=50 Q=250 Q  =  500 

0  1.  1.          1. 

1  0.30373  0.30377-     0.30370 

2  -0.75475  -0.75486  -0.75486 

3  -0.69072  -0.69086  -0.69095 

4  0.23217  0.23227     0.23227 

5  0.63317  0.63357     0.63358 

6  0.15702  0.15717     0.15718 

7  -0.31184  -0.31222  -0.31224 

8  -0.22189  -0.22225  -0.22226 

9  0.05648  0.05660     0.05660 

10  0.11777  0.11807      0.11808 

11  0.02596  0.02604     0.02605 


146 


VI.   CONCLUSION 

Digital  non-recursive  filters  have  been  extensively 
developed  for  the  implementation  of  "finite  impulse  response" 
(FIR)  filters  using  software  on  digital  computers.   In  recent 
years,  new  CTD  (charge  transfer  device)  tapped  delay  lines 
make  it  possible  to  implement  real  time  FIR  filters  using  these 
hardware  devices.   However,  because  the  signal  is  sampled 
analog  instead  of  digital  and  because  the  device  principle 
and  technology  are  new,  careful  investigation  is  needed  for 
the  successful  development  of  these  hardware  sampled  analog 
signal  processors.   In  this  thesis,  the  new  Reticon  TAD-12 
tapped  delay  line  with  external  programmable  tap  weights  is 
investigated  in  detail.   In  addition  to  a  thorough  evaluation 
of  its  device  performance,  the  properties  of  three  sampled 
analog  FIR  filters  using  TAD-12  are  investigated.   They  are 
the  prewhitening,  dewhitening  and  bandpass  filters. 

Both  the  frequency  response  of  TAD-12  devices  have  been 
evaluated  for  sampling  frequency  up  to  2.5  MHz.   It  was  found 
that  the  output  showed  the  following  limitations: 

1.  Frequency  response  limitation: 

Tap  outputs  generally  roll-off  faster  than  predicted 
by  the  effect  of  sample  and  hold.   The  roll-off  depends  on 
the  location  of  the  tap  and  its  tapping  resistance. 

2.  Time  response  limitation: 
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Output  at  one  tap  generally  takes  more  than  the  delay- 
period  between  adjacent  taps  to  settle  down  and  affect  the 
output  of  the  following  taps.   The  settling  time  depends  on 
the  location  of  the  tap,  its  tapping  resistance  and  the 
sampling  frequency. 

3.  Non-uniformities: 

Tap  outputs  are  not  uniform  among  the  taps.   The  non- 
uniformity  varies  with  the  location  of  the  tap,  its  tapping 
resistance  and  signal  frequency. 

4.  Loading  effect: 

Tap  outputs  vary  with  tapping  resistance.   However, 
it  should  be  pointed  out  that  these  limitations  are  large  at 
higher  sampling  frequencies  and  signal  frequencies.   At 
sampling  frequencies  less  than  100  KHz  and  signal  frequencies 
less  than  10  KHz,  TAD-12  behaves  satisfactorily. 

The  properties  of  three  sampled  analog  FIR  filters  using 
TAD-12  have  been  evaluated  following  three  design  procedures. 

1.  Tap  resistances  are  selected  using  the  same  design 
procedure  developed  for  digital  non-recursive  filters. 

2.  Tap  resistances  are  corrected  using  the  frequency  response 

3.  Tap  resistances  are  corrected  using  the  impulse  time 
response. 

It  was  found  that  the  device  limitation  causes  considerable 

deviation  of  the  filter  performance  from  theoretical  prediction 

especially  when  the  sampling  frequency  is  high.   Furthermore, 

the  deviation  is  relatively  smaller  in  bandpass  filter  and 

becomes  worse  in  dewhitening  filter  followed  by  prewhitening 

filter.   This  trend  is  the  result  of  the  relative  magnitudes 
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of  the  tap  outputs  with  respect  to  the  noises  caused  by  the 
device  limitations. 

These  problems  encountered  in  the  Ret  icon  TAD- 12  tapped 
delay  line  just  showed  that  it  is  still  in  its  developing 
stage.   Improved  results  can  be  expected  when  the  causes  of 
the  device  limitations  are  better  understood  and  corrected. 
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APPENDIX    A 


C ******  4**************************************************** 
C       FFCGFA^    OF    CESIGNING    NQNRECURSIVS    PREWHIT5NING    FILTER    BY 
C       UNCCh    METHOD 

Q **************************************************** ******* 
C 

c 

C   THIS  PROGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  ANC  PHASE  OF 

C   FIR  CIGITAL  FILTER  WITH  TRANSFER  FUNCTION 

C   H(Z)  =  C(0)+C(l)*Z**(-l)+C(2)*Z**(-2)+....+C(M )*Z**(-M) 

C   C(M)   FINITE  IMPULSE  RESPONSE 

C   A      SLOPE  CF  FREQUENCY  AMPLITUDE  RESPONSE  FRCM  FO  TO 

C         HIGHEST  FREQUENCY  FP 

C   R      CONSTANT  FREQUENCY  AMPLITUDE  PSEPONSE  FRCM  0  TO  FC 

C   T  =  1.  (NORMALIZED  DELAY  UNIT)  AND  FP=0.5*(1/T)  AND 

C   E  =  FO/FP 

C 

C 

CCNPLEX  ZINV,G,P 

CCNFLEX  CEXP.CMPLX 

DIMENSION  F(150),C(30),W(30),D(20),DD(30),CM20), 
lAfP(15C),PHASE(150) 

UTEGER*4   ITB(12)/12*0/ 

REAL*4     RTB(28)/28*0.0/ 

CATA  R/0.1/,A/l./,E/119.048/,N/5/ 

PI»2. 14159265 

CCNVRT=180./PI 

fr=2*N+l 

kFITE(6,500)  M 
5QC  FOPMATllXt 'LENGTH  OF  FILTER  IS:', 14///) 
C 

C   TC  CCNPUTE  FINITE  IMPULSE  RESPONSE  ANC  WINDOW  FUNCTION 
C 

WCM1)=1. 

C(N+1)=R+A/(2*PI )*(.5-l./(2.*E) )**2 

CC  600  K=1.N 

C(K)=A*(CCS(K*PI >-CCS(K*PI/E)  )/ ( 2.*PI** 3*K**2  ) 

C(N+1-K)=D(K) 

C  <N  +  1+K)  =  C(N+1-K) 

DC(K)=0.54+0.46*CQS(2.*PI*K/M) 

W(N+1-K  )=DD(K) 

MN  +  1+K)  =  W(N  +  1-K) 
600  CCNTINUE 

WFITE(6,650)(C( I) ,I=1,M) 
650  FCFMAT(//2X, 'COEFFICIENTS  ARE : ' /8 ( IX, F10. 5  )  ) 

CNAX=A6S(C(N+1)  ) 

CC  700  1=1, M 

IF(ABS(C(I  D.GT.CMAX)  CMAX=  ABS  (  C  (  I  )  ) 
7CC  CCNTINUE 

CO  750  1  =  1, M 
C 

C   NCRPALIZE  FIR 
C 

C(I)=C(I)/CMAX 
750  CCNTINUE 

WPITE(6,800)(C( I)  ,1=1, M) 
800  FCPMAT(//2X, 'NORMALIZED  COEFFICIENTS  ARE : • /£( IX, F 10.5) ) 

WRITE  <  6,900)  (W(I  ),I=1,M) 
900  FCFMAT(//2X, 'HAMMING  WINDOW  ARE : • /8( IX, FIG. 5 ) ) 

CC  9  1=1. M 
9  CMI)=C(I)*W<I  ) 

W  P ITE  (6,1000  )  (CWU),  I  =  1,M) 
1000  FCFMAT(//2X,' COEFFICIENTS  IMPROVEC  BY  HAMMING  ARE* 
1/8(1X,F10.5)) 

WFITE(6,50) 
50  FCRMAT< ///5X, • FREQUENCY ', 8X ,♦ AM FL ITUDE • , 13X  ,» PHASE • ///  ) 

CC    10    J  =  l,101 

C       F    IS    FREQUENCY    IN    FRACTIONS    OF    NYQUIST    FREQUENCY 

150 


FtJ)=0.G05*FLCAT(J-l) 
ZINV=CE>P(CMPLX(0.,-2.*FC J)*PI )  ) 

PsCMPLXCl.fO.) 

G*CPPLX(0.,0.) 

CC    6    1=1, M 

G=G4P*CW(I  ) 
6    P=P*ZINV 

AfP(J)=CABS(G) 

XsREAUG) 

>=AIMAG(G> 

IF(ABS(X).GT.l.E-8)    GO    TO    60 

PMSE(J)=SIGN(90.,Y) 

GC    TO    8 
60    PPASE(J)=ATAN2(Y,X1*CGNVRT 

8    WPITE<6,20)    F( J), AMP ( J ) , PHASE < J  ) 
20    FCPMAT(5X,F6.3,4X,E16.8,4X,E16.8> 
10    CCNTINUE 

ECLIVALENCE(TITLE,RTB<5)) 

RE4L*8      TITLEU2) 

PEAC(5,300)    TITLE 
300    FCPMATUA8) 

C*LL    CRAWP(101,F,AMP,ITB,RTE) 

PE4C(5,200)    TITLE 

C*LL    DRAHPClOlf F, PHASE, ITB,RTB) 

STCP 

Efvc 

//GC.SVSIN  CC  * 

AMPLITUCE  FREQUENCY  RESPONSE  OF  PREWHITENING 
FILTER  /  LENGTH  =  11  /  HAMMING  WINDOW 
PHASE  RESPONSE  OF  PREWHITENING  FILTER 
LENGTH  *  11  /  HAMMING  WINDOW 
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APPENDIX  6 


C** *************************************************** ****** 

C   FPGGFAM  OF  CESIGNING  CWF  FROM  GIVEN  IMPULSE  RESPONSE  OF 

C   FWF  6Y  FRECUENCY  SAMPLING  METHOD 

C   THE  FRESCRIBEC  CHARACTERISTIC  AND  DESIGNED  FRECUENCY 

C   CHARACTERISTIC  OF  DWF  ARE  PLOTTED  GN  THE  SAME  GRAPH 

C**  ******  ****************************  ************  *********** 

C 
C 

C  TFIS  FROGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  ANC  PHASE  OF 

C  FIR  CIGITAL  FILTER  WITH  TRANSFER  FUNCTION 

C  K2)  -    C(0  )+C(l  )*Z**<-1  )+C(2)*Z**(-2)+....+C(M  )*Z**(-M) 

C  C(M)   FINITE  IMPULSE  RESPONSE 

C 

c 

CCNFLEX    ZINV,G,P,HU50) 

CCNPLEX    CEXPjCMPLX 

ClfENSICN    F{150),C(30),W(30),D(20),DD(30),CM20), 
lAf'P(150),PH'VSE(150),HH(15O) 

INTEGER**      ITB( 12  )/12*0/ 

REAl*4  RTB<28)/28*0.0/ 

FACTOR=I.12 

MM=0 

N»5 

PI*2. 14159265 

CCNVRT=180./PI 

M2*N+1 

WFITE<6,500)  M 
500  FCFMATQX, 'LENGTH  OF  FILTER  IS:1, 14///) 

REAC(5,1)(C(I),I=1,M) 
1  FCFMAT(6F1C5) 
610  WFITE (6,650)(C( I > , 1=1 »M ) 
65C  FGFMAT(//2X, 'COEFFICIENTS  ARE  :  • /8  (IX, F10 .5  )  ) 

C   CCNPLTE  HAMMING  WINDOW 
C 

W(N+1)=1. 
CC  6CO  K  =  1,N 

CC(K)=0.5  4+0.46*C0S(2.*PI*K/M) 
W(N*1-K  )=DD(K) 
W(N+1+K)=W(N+1-K) 
60C  CONTINUE 

CNAX=ABS(C{N+1) ) 
CC    700    1=1, M 

IF(ABS(C(I  H.GT.CMAX)  CMAX=ABS ( C (  I ) ) 
700  CONTINUE 

CC  750  1=1, M 
C 

C   NCPNALIZE  FIR 
C 

C(I)=C(I)/CMAX 
750  CCNTINUE 

V»FITE(6,800)(C(I),I=1,M) 
800  FCPMATC//2X, 'NORMALIZED  COEFFICIENTS  ARE  :  ■ /8( IX, F10.5)  J 
CC  S  I=1.M 
«  CMI)»C(I)*W(I) 

WPITE(6,1000)(CW(I) ,I=1,M) 
1000  FCPNAT(//2X, 'COEFFICIENTS  IMPROVED  BY  HAMMING  ARE' 
1/8<1X,F10.5») 
WFITE(6,50) 
5C  FC F MAT (/// 5X, • F REQU EN CY',8X,« AMPLITUDE' , 13X ,« PHASE •///  ) 
DC  10  J=l,101 
C 
C   F  IS  FREQUENCY  IN  FRACTIONS  OF  NYQUIST  FREQUENCY 


C 


F(J)=0.C05*FL0AT( J-1J 
ZINV=CEXP(CMPLX(0.,-2.*F(J)*PI ) ) 
F=CNPLX(1.,0.) 
G  =  CfPLM0.,0.) 
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CC  6  1=1, M 

IF(MM.SC.l)  GO  TO  11 
G=G*P*CW(I ) 
GC  TO  6 
11  G=G4P*CW(I)*FACTOR 
i    F=F*ZINV 

*MF(J)=CABS(G) 

c 

C   FPECLENCY  SAMPLE  POINTS  CF  DWF 
C 

F(J)=1./G 

HH(J)=CABS(H(J) ) 

X=PEAL(G) 

Y=AIMAG(G) 

IF<ABS(X).GT.l.E-8)  GO  TO  60 

PFASE(J)=SIGN(90.,Y) 

GC  TO  8 
6C  PHASE(J)=ATAN2(Y,X)*C0NVRT 

8  WPITE<6,20)  F( J),AMP( J) ,PHASE( J) 
2C  FCFMAT(5X,»:6.3,4X,E16.8,4X,E16.8) 
10  CCNTINUE 

ECUIVALENCE(TITLE,RTB(5) ) 

RE*L*8   TITLE(12) 

IF(MM.EC.l)  GO  TO  301 

GC  TO  302 

301  ITE(1)=2 
ITE(2)=C 

C4LL  CRAWP(101,F,AMP, ITB,RTB) 

302  L=N  +  1 
C 

C   CF  CWF 


C 

CC  1500  1=1, L 

K  =  I-1 

GG  =  CMPLMO.,0.) 

DC  1100  J=l,101 

Z=CEXP(CMPLX(0.,2.*PI*F(J)*K) ) 

GG=GG+Z*H< J) 
11CC  CCNTINUE 

C(I)=GG/101. 

C (N+l-I  )=C(I) 
1500  CCNTINUE 

N  IW  —  N  M  +  1 

IF(MM.EC.2)  GO  TO  2100 

PE*C<5,200)    TITL5 
3CC    FCPMTUA8) 

ITE(1)=1 

ITE(2)=1 
C 

C   FLCT  THE  PPESCRIBEO  FREQUENCY  CHARACTERISTIC  CF  CWF 
C 

C   LSING  INVEPSE  DISCRETE  FCURIER  TRANSFORM  TO  CALCULATE  THE 
C   IMPULSE  RESPONSE  OF  CWF 

CALL  CRAWP(101,F,HH,ITB,RTB) 
C 

C   GC  BACK  TO  PLOT  THE  AMPLITUDE  FREQUENCY  RESPONSE  OF  DWF 
C 

GC  TO  610 
2100  STCF 

ENC 
//GC.SYSIN  CO  * 

-C. 00091    0.0000    -0.01219    0.0000    -0.21427    1.0000 
-0.21427    CO       -0.01219    0.0       -0.00091 
AMPLITUDE  FRECUENCY  RESPONSE  OF  DEWHITEMNG 
FILTER  /  LENGTH  =  11  /  HAMMING  WINDOW 
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APPENDIX    C 


C*******& ******************************* ******************** 

C*    FPCGF/S  M    OF    CESIGNING    NONFECURSIVE    BANDPASS     FILTER 

C ***************************************** ****************** 

C 

c 

C  THIS  FROGRAM  CALCULATES  ANC  PLOTS  AMPLITUDE  ANC  PHASE  OF 

C  FIR  CIGITAL  FILTER  WITH  TRANSFER  FUNCTION 

C  F(2)    =    C(0)+C<l)*Z**(-l)+C<2)*Z**(-2)+....+C(M  )*Z**(-M) 

C  C(M>       FINITE    IMPULSE    RESPONSE 

C  A         AMPLITUDE    OF    BANDPASS 

C  Fl    €    F2         I    CRITICAL    FREQUENCIES 


CCfrFLEX    ZINV,G,P 
CCMFLEX    CEXP.CMPLX 

DIMENSION    F(150),C<30),W(30),D<2G),DD<30)  ,CM3C)  , 
^F<150  ),PHASE(150) 
INTEGER*4      ITB(  121/12*0/ 
REAL*4  RTB(28)/28*0.G/ 

CATA    A/l./,Fl/0.1980/,F2/0.2020/,N/5/ 
PI=3. 14159265 
CCNVRT-180./PI 
M2*N+1 

WPITE<6,500)    M 
50C    FCFMAT(1X, 'LENGTH    CF    FILTER    IS:', 14///) 


w  <n+i)=i. 

C(N+1)=2.*A*(F2-F1) 

CC  600  K=1,N 

C  <K)=A*(SIN(2.*PI*K*F2)-SIN(2.*PI*K*F1)  )/<K*PI) 

C(M1-K)  =  D(K) 

C  <N+1-H<)  =  C<N+1-K) 

CC(K)=0.54+0.46*COS(2.*PI*K/M) 

W<N+1-K)=CC(K) 

W(N*1+K  )=W(N+1-K) 
600  CCNTINUE 

WPITE(6,650XC( I),I=1,M) 
650  FCPMAT(//2X, 'COEFFICIENTS  ARE: • /fi ( IX, F10. 5  )  ) 

CNAX=ABS(C(N+1) ) 

DC  700  I  =  1?M 

IF(ABS(CU  M.GT.CMAX)  CMAX=ABS ( C (  I  )) 
700  CCNTINUE 

CC  750  1=1, M 

C      NORMALIZE    FIR 
C 

C  <I)=C(I)/CMAX 
7SC    CCNTINLE 

WFITE(6,80C)(C( I) ,1=1, M) 
80C    FCFMAT(//*2X, 'NORMALIZED    COEFFICIENTS    ARE:  •  /€(  IX,  F10.5)  \ 

WRITE    (6,900)    <W<I),I=1,M) 
900    FCPNAT(//2X, 'HAMMING    WINDOW    ARE  :  •  /8 ( IX, *=10. 5 )  ) 
DC    9    1=1, M 

?  CMi)=cm*wm 

WPITEt6,1000)(CW(I),I=l,M) 
1000  FCFMAT(//2X, 'COEFFICIENTS  IMPROVED  BY  HAMMING  ARE' 
1/6(1X,F10.5)> 
WPITE(6,50) 
50  FORMAT (///5X,'FREQUENCY',8X,« AMPLITUDE'  , 13X  , ' PHAS E  • ///  1 
CC  10  J  =  l,  101 

C   F  IS  FREQUENCY  IN  FRACTI  CNS  OF  NYQL'IST  FREQUENCY 


C 


F(J)=0.01*FL0AT(J-1) 
ZIN'V=C£XP<CMPLX(0.,-2.*f=(  J)*PI  )  ) 
F=CMPLX(1.,0.) 
G=CNFLX(0.,0.) 
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cc 

Q  =  G 
F*F 
ANP 
X  =  F 
V^ 
IF( 
PH/ 
GC 
PH* 
WPI 
FGF 
CCN 
ECU 
RE* 
FE* 
FCF 
C4L 
RE* 
CAL 
STO 
ENC 
//GC.SY5I 
AMPLITUDE 
FILTER  / 
PHASE  PES 
LENGTH    = 


6G 

6 

20 

10 


200 


6    1=1, M 

4P*CW(I) 

*ZINV 

(J)=CABS(G) 

EAL(G) 

IMAG(G) 

ABS(X).GT.l.E-S) 

SE(J)=SIGN(90.,Y) 

to  a 


GO    TO    60 


TO    8 

iSE(J)=ATAN2( Y,X)*CONVRT 
TE(6,20)    F(J),AI*P(  J),PHASE(J) 
!^AT(5X, F6.3,4X,E16.8,4X,E16.8 
TINUE 

!IVALENCE(TITLE,RTB(5)) 
L*8      TITLE(12) 
C(5,20G>    TITLE 
*AT(6A8) 

L    DRAWP(101,F,AMP,ITB,RTB) 
C<5,300)    TITLE 
L    DRAWP(101,F, PHASE, ITB,RT 
P 


RTB) 


FREQUENCY  RESPONSE  OF  BANDPASS 
LENGTH  =  11  /  HAMMING  WINDOW 
FCNSE  OF  BAN0PA5S  FILTER 
11  /  HAMMING  WINDOW 
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APPENDIX  C 


C   THIS  FROGRAM  CALCULATES  IMPULSE  RESPONSE  GIVEN  AMPLITLDE 
C   FPEQLENCY  RESPONSE  OF  FILTER 

CCMFLEX  CEXP.CMPLX 

C ^ENSIGN  F(150),C(30),W(30),CW(2C),DO(30),I-(150) 

N  =  5 

P 1=3. 14155265 

M2*N+1 

WFITE(6,5QC)  M 
5CC  FCPMATC1X, 'LENGTH  OF  FILTER  IS:1, 14//) 
C 

C   CCMPLTE  HAMMING  WINDOW 
C 

MN  +  1)=1. 

DC  600  K=1,N 

CC(K)=0.54+0.46*C0S(2.*PI*K/M) 

W(N+1-K)=DD(K) 

W(N-H+K)  =  W(N  +  1-K) 
600  CCNTINUE 

PEAC(5,7C0>  (H(J),J=1,101) 
7CC  FCFMAT<6F10.5) 

WRITE<6,750)  ( H( J ) , J  =  1, 10 1 ) 
75C  FGFMATC//2X, 'DATA  ARE  :  •  ,/// (8 ( 1 X  ,  F10.5 I,/ ) ) 
C 

C   NCFMLIZE  MJ) 
C 

FNAX=H<89) 

CC  2000  J=l,101 

IF(H(J).GT.HMAX)  HMAX=H(J) 
2000  CCNTINUE 

CC  2100  J»l»101 

H(J)=H( JJ/HMAX 
2100  CCNTINUE 
C 

C   LSING  INVERSE  DISCRETE  FCURIER  TRANSFORM 
C   TC  CALCULATE  THE  IMPULSE  RESPONSE 
C 

L=M 

CC  1500  1  =  1, L 

K  =  I-1 

GG=CMPLX(Q.,0.) 

DC  1100  J  =  l,101 

F(J)=0.C1*FLQAT(J-1) 

Z=CEXP(CMPLX(0.,2.*PI*F(J)*K) ) 

GG=GG+Z*H(J) 
HOC  CCNTINUE 

CCI)=GG/101. 
1500  CCNTINUE 

V»PITE(6,650)(C(I)  ,1=1, M) 
650  FCRMAT(//2X, 'COEFFICIENTS  ARE: ■ /8 ( IX, F10. 5 )  ) 

C   NCFMALIZE  FIR- 
C 

CMX=ABS(C(N+1) ) 

CC  800  1  =  1, M 

IF(ABS(C(I>).GT.CMAX)  CMAX=ABS( C ( I ) ) 
800  CCNTINUE 

DC  $00  I  =  1,M 

c(ii=ci  d/cmax 

90C    CCNTINUE 

WRITE(6,85GMC( I) ,1=1, M) 
850    FCFMAT(//2X, 'NORMALIZED    COEFFICIENTS    ARE  :  • /£(  IX, F1Q.5) 
CC    9    I  =  1.M. 
c    CK(I}»Cm*W(I) 

WRITE  (6,1000  MCW< I)  ,1  =  1, M) 
1000    FCRMAT( //2X, 'COEFFICIENTS    IMPRCVEC    BY    HAMMING' 
l/8(  lXrF10.5) ) 
STCF 
ENC 
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APPENDIX  E 


C**** ********* *********************************** *********** 
C   FPCGFAM  FOR  PLOTTING  100  CALCULATED  AND  MEASURED  POINTS 

C   CF  ANY  NCNFECURSIVE  FILTER  WITH  GIVEN  IMPULSE  RESPONSE 
£******** *************************** ************************ 

C 

c 

C   THIS  FPCGRAM  CALCULATES  AND  PLOTS  AMPLITUDE  AND  PHASE  OF 

C   FIR  CIGITAL  FILTER  WITH  TRANSFER  FUNCTION 

C   MZ)  =  C<0)+C(l)*Z**(-l)  +  C(2)*Z**<-2)+....+C(M)*Z**<-y) 

C   C(N)   FINITE  IMPULSE  RESPONSE 

C 

C 

CCNFLEX    ZINV,G,P 
CCNFLEX    CEXP.CMPLX 

DIMENSION    F(150),C(30),AMP( 150)  ,  PHASE (150) ,B(150)  , 
IE E < 150 ),B A (150), AM  PL (150),  AMPLTU50) 
1AVFL(150) ,AMPLT(15C) 
INTEGERS         ITB ( 12  )/12*C/ 
PEAL*4  RTB(28)/28*0.0/ 

PI=2. 14159265 
N  =  ll 

CCNVRT=180./PI 
PEAC(5,5C0)(C( I),I=1,M) 
50C    FCRNAT(8F10.5) 

RE4C(5,5O0)    (B(K)  ,K=1,101) 
fcPITE(6,850)    (B(K),K=1,101) 
85C    FCFMAT(//2X,»DATA    ARE : • ,/// (8 ( IX , F10.5 ) , / )  ) 
BNAX=B( 69) 
DC    1000    J=l,101 
10C0    I F ( e( J) .GT.BMAX)     BMAX=B(J) 
DC    HOC    J=l,101 
EE(J)=B(J)/BMAX 
HOC    EMJ)=20.*ALOG10(B8(J)  ) 
WFITE(6,8GC)(C( I),I=1,M) 
80C    FCPMAT(//2X, 'NORMALIZED    COEFFICIENTS    ARE  : '/£( IX, F 10.5)  ) 
WPITE(6,50) 
50    FCFMAT(///5X, • FREQUENCY1, 8X , • AM PL ITUDE ■ , 13 X ,' PHASE •/// ) 
CC    10    J=l,101 

c 

C   F  IS  FREQUENCY  IN  FRACTIONS  OF  NYQUIST  FREQUENCY 


C 


F(J)=Q.C05*FL0AT(J-1) 

2INV=CEXP(CMPLX(0.,-2.*F( J)*PI ) ) 

F=CNPLX(1.,0.) 

G=CNPLX(0.,0.) 

CC  6  1=1, M 

G=G+P*C(I) 
6  F=P*ZINV 

A^F(J)=CABS(G) 

>=PEAL(G) 

Y=AIMAG(G) 

IF(ABS(X).GT.l.E-8)  GO  TO  60 

PFASE(J)=SIGN(90.,Y) 

GC  TO  8 
60  PFASE(J)=ATAN2(Y,X)*CCNVRT 

6  WPITE(6,20)  F( J) ,AMP( J) ,PHASE(J) ,E(J) 
20  FCFMAT(5X,F6.3,4X,E16.8,4X,E16.8,6X,E16.8) 
10  CCNTINUE 

ANAX=AMF(51) 

CC    2C00    JaltlOl 

IF(AMP( J).GT.AMAX)    AMAX=AMP(J) 
20C0    CONTINUE 

CC    2100    J*1,101 

A^PL( J)=AMP( J)/AMAX 
2100    ANPLT(J)=20.*AL0G10(AMPL( J) ) 

iTe  ci >  =  i 

I7E(2)=0 
ECUIVALENCE(TITLE,RT3(5)) 
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REAl*8         TITLEQ2) 

RE£C(5,7)    TITLE 

FCFN*T(6*8) 

C*LL    CRAWPQ01,F,AMPLT,ITB,RTB) 

ITB(1)=3 

ITE(2)=3 

C4LL    CRAWP(101tF,BA»ITB,RTB) 

STCP 

ENC 
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APPENDIX  F 


TAPPING  RESISTANCES  FOR  NON-RECURSIVE  FILTER 


tap 


■k 


I 


out 


FIGURE  P.1     Coefficient  Implementation 


The   contribution   of   a   single   tap,V        to    the   output,   V 
is    the  product   of  V,    and   the   coefficient,   A,  ,    or 


out 


Rf 


r  +  Ri 


The  expression  indicates  that  R^  is  inversely  proportional 
to  Afc,  therefore  the  largest  coefficient  will  correspond  to 

the  minimum  resistance,  R  •„,  so 

'   mm 


Rf 


max 


r  +  R 


mm 


Solving  this  equation  for  R^  and  substituting  that  result 
into  the  A,  expression  yields: 


A, 


r  +  R, 
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or, 

„    Amax 

Rk  =  ~K]T   ^  Rmm)  -  r 

For  the  TAD-12,  the  output  impedance  for  a  tap,  r,  is  about 
5  Kilo-ohms. 
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